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The  primary  goal  of  this  contract  is  to  pursue  research  studies  of 
tunable  optical  sources,  to  develop  coherent  sources  in  new  spectral  re- 
gions, and  to  conceive  and  to  investigate  new  nonlinear  optical  devices 
of  interest  to  the  Air  Force.  During  this  contract  our  efforts  have  been 
concentrated  on  two  primary  projects:  experimental  studies  of  possible 

metal  vapor  excimer  laser  systems,  and  a unique  cesium  vapor  infrared  to 
visible  frequency  converter  for  imaging. 

Early  in  this  program  a number  of  metal  vapor  excimer  laser  systems 
were  investigated  experimentally  with  the  hope  of  demonstrating  net  opti- 
cal gain.  Based  on  a number  of  factors,  it  was  hoped  that  such  excimer 
systems  might  prove  useful  as  high  efficiency,  high  power  density  ampli- 
fiers or  lasers  in  the  visible  or  near  infrared  portion  of  the  spectrum. 
Because  of  the  lack  of  detailed  knowledge  of  the  interatomic  potential  of 
such  systems,  our  investigation  was  primarily  experimental.  Na-Ze,  Li-Xe, 
and  Na-Na  systems  were  studied,  but  no  net  gain  was  observed  in  any  case. 

In  addition,  these  experiments  revealed  a number  of  severe  technical  prob- 
lems related  to  the  construction  of  high  temperature,  high  pressure  cells 

containing  reactive  metals.  In  addition,  the  uniform  optical  pumping  of 

-/ 

» 

such  highly  absorbing  and  dispersive  atomic  vapors  proved  very  difficult. 

In  view  of  these  problems,  along  with  the  very  small  optical  gains  availabl 


we  subsequently  devoted  our  entire  efforts  to  the  cesium  infrared  image 
up-converter  project. 

The  IR  image  up-conversion  technique  employs  a three-photon  sum  pro- 
cess in  cesium  vapor.  High  efficiencies  are  achieved  by  two-photon  pump- 
ing a non-allowed  transition  in  order  to  produce  a resonantly  enhanced 
nonlinearity;  using  this  large  nonlinearity,  image  up-converters  having 
high  efficiency  and  resolution  can  be  constructed  using  only  moderate  pump 

powers.  The  solid  state  pumping  laser,  Nd: lanthanum  beryllate,  has  a nat- 

2 2 

ural  two-photon  coincidence  with  the  cesium  6s  S-7s  S transition,  resulting 
in  a simple,  practical  system  with  a number  of  potential  applications  for 
communications,  reconnaissance,  and  the  imaging  of  rapid  events. 

An  evaluation  of  IR  imaging  technologies  sponsored  by  the  Air  Force 
* 

Avionics  Laboratory  concludes  that  the  overall  performance  of  such  metal 
vapor  techniques  is  considerably  superior  to  crystal  systems,  and  is  com- 
parable to  proposed  direct  IR  detection  methods.  A particular  advantage 
of  the  nonlinear  technique  is  the  high  instantaneous  conversion  efficiency 
achieved;  this  makes  time  resolved  imaging  of  fast  events  possible,  a ca- 
pability notably  lacking  in  direct  detection  devices  which  require  long 
integration  times.  In  addition,  the  relative  timing  of  an  illuminating 
source  and  the  up-converter  pump  can  be  used  for  range-gated  imaging,  both 
to  determine  range,  and  to  reduce  interference  from  undesired  objects. 

Other  attractive  features  of  these  systems  include  scalability  to  large 

* 

Kamala  S.  Krishnan  and  John  S.  Ostrem,  "Evaluation  of  Technologies 
for  Infrared  Imaging,"  Stanford  Research  Institute  Project  ISE  U0U7  for 
AFAL/TEO  Contract  F33615-75-C-IIU2  (January  1976). 


apertures  and  collection  angles,  high  transmission  throughout  the  IR,  and 
the  ability  to  withstand  high  incident  power  densities  without  damage. 

Our  experimental  work  has  demonstrated  a number  of  these  device  char- 
acteristics. We  have  recently  demonstrated  the  up-conversion  of  . . 9 p IR 
images  to  U550  A in  cesium  vapor.  A power  conversion  efficiency  of  20fo 
with  over  1000  resolvable  spots  was  achieved  using  a pump  power  of  only 
8 kW . In  addition,  an  extensive  theoretical  study  of  the  factors  influ- 
encing image  up-converter  efficiency,  resolution,  and  image  aberrations 
has  been  performed.  This  work,  including  details  of  the  experimental  mea- 
surements, and  optimized  designs  for  higher  performance  up-converters,  is 
included  as  Appendix  A. 

The  publications  resultihg  from  work  under  this  contract  are  listed 
in  Appendix  B;  previous  related  contracts  and  their  publications  are  listed 
in  Appendix  C.  Portions  of  this  work  were  jointly  supported  by  the  Army 
Research  Office  and  the  Joint  Services  Electronics  Program. 
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ABSTRACT 


Two-photon-resonant  frequency  converters  are  shown  to  offer  an 
attractive  method  for  up-converting  infrared  images  into  the  visible 
and  ultraviolet.  It  is  found  that  for  this  application  they  are  far 
superior  to  crystal  up-converters.  Experimental  results,  including 
the  first  images  produced  with  this  technique,  are  presented. 

A theoretical  analysis  of  two-photon-resonant  processes  is  per- 
formed, both  for  the  cases  of  Q-switched  and  mod  —locked  lasers.  Gen- 
eral results  are  obtained  for  efficiency,  two-photon  absorption  proba- 
bilities, and  the  Kerr  contribution  to  the  index  of  refraction.  A 
new  method  for  measuring  linewidths  of  non-allowed  transitions,  based 
on  the  results  of  this  analysis,  and  experimental  results  demonstrating 
this  method,  are  presented. 

Finally,  harmonic  generation  and  mixing  processes  at  high  in- 
tensities are  discussed.  It  is  shown  how  the  resonances  of  nonlinear 
processes  are  intensity-dependent.  The  effect  of  these  "frequency 
shifts"  on  the  efficiency  of  VUV  and  soft  x-ray  generators  is  dis- 
cussed . 
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CHAPTER  I 


INTRODUCTION 


This  dissertation  presents  a theoretical  and  experimental  study 
of  the  feasibility  of  using  two-photon-resonant  frequency  conversion 
processes  in  metal  vapors  as  a technique  for  up-converting  infrared 
images  into  the  visible  or  near  UV.  The  experimental  results  include 
(the  first)  images  obtained  using  this  technique.  The  superiority  of 
the  technique  over  up-conversion  in  nonlinear  crystals  is  clearly 
demonstrated . 

Also  presented  is  a theoretical  analysis  of  the  dependence  of  the 
efficiency  of  two-photon-resonant  frequency  converters  on  the  laser 
and  atomic  linewidths.  A new  technique  for  measuring  linewidths  of 
non-allowed  transitions,  based  on  the  results  of  this  analysis,  and 
which  has  been  experimentally  demonstrated,  is  also  discussed. 

Finally,  a theoretical  analysis  is  presented  of  the  interaction 
of  atoms  with  very  intense  laser  fields  as  are  used  for  the  generation 
of  wavelengths  in  the  VUV  range  of  the  spectrum. 

The  concept  of  up-conversion  in  nonlinear  crystals  is  a number  of 

1.2 

years  old.  ’ The  merits  of  such  systems  for  infrared  imaging  applica- 
tions were  discussed  in  a review  paper  by  F.  Milton.  Two  severe 
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drawbacks  of  such  systems  are  the  very  high  powers  required,  together 
with  the  poor  optical  and  mechanical  properties  of  most  nonlinear 
crystals.  Recently,  however,  a new  technique  for  frequency  up-con- 

, k-6 

version  has  been  demonstrated,  which  makes  use  of  the  strong  non- 
linearities  of  alkali  metal  vapors  when  the  pump  radiation  is  tuned 
to  a two-photon  resonance.  Because  of  the  resonant  nature  of  the 
interaction,  the  powers  required  for  high  efficiency  high  resolution 
imaging  are  typically  more  than  two  orders  of  magnitude  less  than 
those  required  when  making  use  of  the  best  available  crystals,  while 
damage  problems  do  not  exist.  The  price  to  pay  for  these  improvements 
is  that  a pump  laser  is  required  which  either  has  an  accidental  reso- 
nance with  a non-allowed  transition  in  an  alkali  metal,  or  which  can 
be  tuned  to  such  a resonance.  Also,  the  choice  of  IR  wavelengths  that 
can  be  up-converted  is  limited  because  the  generated  sum  frequency  has 
to  come  fairly  close  to  an  allowed  transition  of  the  vapor.  In  Chap- 
ter III,  infrared  imaging  with  such  an  up-converter  is  discussed  in 
detail.  The  various  aspects  involved  are  explained  using  the  example 
of  3 p.  up-conversion  in  atomic  cesium,  for  which  a very  convenient 
pump  source  is  available. 

The  efficiency  of  a two-photon-resonant  frequency  converter  is, 
in  the  most  general  case,  a function  of  the  laser  spectrum,  the  pres- 
sure broadened  linewidth  and  the  Doppler  width.  General  expressions 
are  derived  for  this  efficiency  for  two  cases.  In  the  first  case  the 
modes  of  the  laser  are  assumed  to  be  independent  (Q-switched  laser  , 
while  in  the  second  case  they  are  assumed  to  be  all  zero  (mode-locked 
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laser) . 


The  results  of  the  first  case  suggest  a new  simple  way  for 


ft 


measuring  pressure  broadened  linewidths,  requiring  only  a tunable  laser 
that  may  be  several  orders  of  magnitude  wider  than  the  linewidth  being 
measured.  This  has  to  be  contrasted  with  other  existing  techniques 
which  all  require  a laser  narrower  than  the  atomic  linewidth.  ‘ The 
technique  also  is  able  to  measure  linewidths  which  are  more  than  an 
order  of  magnitude  narrower  than  the  Doppler  width.  The  analysis  of 
two-photon-resonant  processes,  together  with  experimental  results  for 
the  pressure  broadening  by  He  of  the  Js  S-3d  D transition  of  Na,  are 
presented  in  Chapter  II. 

Chapter  IV,  which  is  unrelated  to  the  rest  of  this  dissertation, 
discusses  harmonic  generation  mixing  processes  which  require  very 
high  intensities.  These  processes  offer  a promising  technique  for  the 
generation  of  radiation  in  the  wavelength  range  of  a few  hundred  ang- 
stroms . At  the  very  high  intensities  involved,  the  resonances  of  the 
nonlinear  process  are  not  the  same  as  in  the  weak  field  case,  but  thev 
become  dependent  on  the  power  of  the  laser  s . Although  these  "fre- 
quency shifts"  can  in  principal  be  compensated  for,  e.g.,  by  tuning  to 
the  new  resonance  frequency  in  the  case  of  a two-photon-resonant  fre- 
quency converter,  this  is  often  very  hard  in  practice  because  of  the 
pulse-to-pulse  fluctuations  of  the  lasers  used  for  such  experiments. 
At  the  present  time  an  effort  is  underway  in  the  laboratory  of  Dr. 
Harris  to  generate  625  X radiation  in  He  by  a third  order  mixing  pro- 
cess of  the  seventh  harmonic,  the  ninth  harmonic,  and  the  fundamental 
of  the  1 .06U  p line  of  the  Nd:YAG  laser. 
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CHAPTER  II 


LINEWIDTH  DEPENDENCE  OF  TWO -PHOTON -RESONANT  PROCESSES. 
APPLICATION  TO  THE  MEASUREMENT  OF  LINEWIDTHS  OF 
ELECTRIC  DIPOLE  FORBIDDEN  TRANSITIONS 

A.  INTRODUCTION 

Two-photoa-resonant  processes  have  been  shown  to  offer  a very  ef- 

L — t 

ficient  technique  for  the  generation  of  UV  as  well  as  IR  radiation. 
They  also  offer  a way  for  populating  atomic  or  molecular  energy  levels 
via  two-photon  absorption  to  some  higher  level,  followed  by  rapid  de- 
cay to  the  level  of  interest. ^ In  the  first  application  these  pro- 
cesses are  attractive  because  they  suffer  far  less  from  absorption  and 
dispersion  at  the  pump  and  the  generated  wavelength  than  do  single-  or 
three-photon-resonant  processes.  In  the  second  application  they  are 
sometimes  preferred  over  direct  population  of  the  pertinent  level  be- 
cause of  the  possibility  to  excite  a certain  small  volume  inside  the 
medium  by  appropriate  focusing  of  the  pumping  radiation,  and  also  be- 
cause the  longer  wavelengths  required  are  often  more  easily  available. 
In  this  chapter  the  calculation  of  such  processes  is  consider.!  in  de- 
tail. First  the  basic  relations  are  presented  from  which  the  effi- 
ciency of  frequency  converters,  two-photon  absorption  probabilities, 
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and  the  quadratic  Kerr  effect  can  be  calculated.  Obviously  there  is  a 
large  number  of  cases  that  can  be  considered,  depending  on  the  proper- 
ties of  the  incident  radiation,  i.e.,  continuous  or  discrete  spectrum, 
independent  or  coupled  modes,  etc.  Two  cases  are  treated  in  detail. 

In  the  first  one,  the  incident  radiation  is  assumed  to  consist  of  a 
set  of  discrete  modes  with  independent  phases.  The  second  case  is 
that  of  a time-bandwidth  limited  pulse  with  a continuous  spectrum. 
Sources  which  produce  pulses  having  a spectrum  of  the  first  kind 
usually  have  pulse  lengths  that  are  long  compared  to  the  atomic  de- 
phasing times,  and  as  a consequence  build-up  times  can  and  will  be 
neglected  in  these  calculations.  For  the  second  case,  which  is  typ- 
ically that  of  a very  short  pulse  switched  out  of  a mode-locked  train, 
an  exact  solution  is  presented,  taking  into  account  the  transitory  be- 
havior. Finally,  based  on  the  results  of  the  first  case,  a new  method 
is  presented  for  measuring  linewidths  of  electric  dipole  forbidden 
transitions,  requiring  for  the  measurement  a laser  source  which  can 
be  several  orders  of  magnitude  wider  than  the  pressure  or  Doppler 
broadened  atomic  linewidth  being  measured. 

B.  BASIC  RELATIONS  FOR  THE  CALCULATION  OF  TWO-PHOTON -RESONANT 
PROCESSES 

1 . The  Nonlinear  Polarization  in  Terms  of  the  Incident  Fields 

The  general  expression  for  the  J>-d  order  nonlinear  susceptibility 

12 

is  given  by  Butcher.  For  the  case  of  two-photon-resonant  processes, 
the  polarization  at  the  third  harmonic  and  the  back-polarization  at  the 
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fundamental  frequency  can  be  calculated  from  the  following  set  of  ex 
pressions : 


where 


X 
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In  the  case  of  J-d  harmonic  generation,  the  contribution  of  the 
last  two  terms  in  the  susceptibility  is  small,  and  we  will  neglect  it 
in  what  follows.  Assuming  that  the  process  is  only  two -photon-reso- 
nant, combination  of  the  expressions  under  (2.1)  gives  the  following 
result : 


P'°!(t;3«>0)  = 


^01^12^23^0 
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+“2)t 


(“*20  +“l  +CU2  + 


where  a summation  over  paths  is  implicit.  Note  that  although  all  the 
calculations  will  be  done  for  third  harmonic  generation,  the  results 
will  be  independent  of  the  field  used  in  the  third  step. 

In  the  case  of  the  back-polarization  at  the  fundamental  frequency, 
the  third  and  fourth  terms  in  the  susceptibility  are  equal  respectively 
to  the  first  and  second  terms.  Assuming  again  that  the  process  is  onlv 
two -photon -re sonant,  the  expression  for  the  polarization  can  be  written 
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where  again  a summation  over  paths  is  implicit.  In  this  case  the 
dominant  term  in  the  summation  over  the  }-d  level  is  the  one  where 
the  3-d  level  is  the  same  as  the  first  level.  In  what  follows  we 
will  only  keep  that  term. 

. . The  Generated  Electric  Field  in  Terms  of  the  Nonlinear  Polarization 

electric  field  is 


The  relation  between  polarization  and  radiates 
given  by  the  wave  equation.  Neglecting  loss: 


1 d E d'“P 
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If  the  polarization  is  expanded  in  a set  of  plane  waves,  a traveling 
wave  equation  can  be  written  for  the  complex  envelope  quantities  of 
each  wave.  In  one  dimension: 


iou  n P 
m m n 


mode  m 


(2.5) 


Because  of  the  linearity  of  the  wave  equation,  the  total  field  is 
found  by  summing  the  contributions  of  the  individual  modes: 


E(z,  t) 
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3 . The  Two-Photon  Absorption  Probability  in  Terms  of  the 
Back -Polarization 


From  energy  conservation: 


- PE  = .NW 


where  P is  the  back-polarization,  N is  the  number  density,  and 

(2) 

W'  is  the  two-photon  absorption  probability  per  second.  From  this 
expression  an  average  absorption  rate  can  be  calculated: 


= 


X lira  — 
T->»  x 


PE  dt 
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Often  it  it:  more  meaningful  to  talk  about  the  total  absorption  proba- 
bility during  a pulse: 
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4 . The  Quadratic  Kerr  Effect  in  Terms  of  the  Back-Polarization 

The  in-phase  component  of  the  back-polarization  gives  a contribu- 
tion to  the  index  of  refraction.  For  infinite  plane  waves,  this  in- 
phase  component  can  be  written  as: 


P(3) 

in  phase 


e0*HE 


(2.10) 


where  E is  the  amplitude  of  the  wave  at  frequency  '»  . The  index 

of  refraction  is  then  given  by  (from  Maxwell's  equations): 


(2.11) 


For  short  pulses,  the  effect  on  pulse  propagation  through  the  medium 
is  found  by  solving  the  T.W.E.: 
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C.  EFFICIENCY  OF  TWO-PHOTON -RESONANT  FREQUENCY  CONVERTERS 

1 . Electric  Field  With  a Discrete  Spectrum  of  Independent  Modes 

We  calculate  the  tripling  efficiency  for  an  incident  electric 
field  of  the  form: 


E ( t ' 


+ ioYi.>)t  - cp  ] 

ra 


(2.13) 


where  the  phases  cp^  are  independent  random  variables.  First  we  will 
calculate  the  polarization  at  the  3_d  harmonic  frequency  using  (2.2'  . 
From  this  we  can  find  the  generated  electric  field  using  the  T.W.E. 
(2.5).  The  efficiency  is  then  calculated  as  the  ratio  of  the  power 
in  the  third  harmonic  field  to  the  power  in  the  fundamental  field. 

Substituting  the  expression  for  the  electric  field  (2.13)  in  the 
expression  for  the  transform  of  the  field  (2.1),  we  find: 
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Substituting  this  in  the  expression  for  the  polarization  (2.2): 
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To  incorporate  Doppler  broadening  in  the  theory  we  note  that,  for 
an  observer  traveling  with  the  incident  electric  field,  atoms  moving 
at  different  velocities  have  different  apparent  resonance  frequencies 
cDq2 (x)  = wqo  + x • F°r  a Maxwellian  velocity  distribution,  the  fre- 
quency deviation  x has  a Gaussian  probability  density  function: 


p(x)  = 


where 
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(2.16) 


There  are  N x p(x)  X dx  atoms  with  (apparent)  resonance  frequencies 
between  “20  + x and  “20  + X + dx  . The  total  polarization  is  found 
by  summing  over  all  velocity  classes.  The  result  is: 
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functions  F (k)  and  F_(k)  are  the  real  and  imaginary  parts  of 
K 1 
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the  plasma  dispersion  function.  Substituting  the  expression  for 
E . (t  in  (2. IT),  the  polarization  can  finally  be  written  as: 
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(2.1  ) 

Using  the  T.W.E.  (.  .3  , the  electric  field  generated  by  each  term 
in  the  sum  (2.1  can  now  be  calculated.  The  total  electric  field  is 
found  by  summing  all  these  contributions.  Assuming  the  same  k-vector 
mismatch  for  all  the  terms  in  the  sum  which  may  not  be  exact  at  in- 
tensities where  the  Kerr  effect  becomes  important,  see  Section  II-E\ 
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the  T.W.E.  for  ail  arbitrary  term  of  the  sum  is: 
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Because  the  phases  (<+  + cp  + Cf  ) for  different  terms  of  the  sum  are 
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independent,  the  time  average  value  of  the  square  of  the  electric  field 
waveform  i;  simply  given  by  one-half  the  sum  of  the  squares  of  the  am- 
plitudes of  the  terms: 
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(Note  that,  since  terms  with  different  p values  have  different  fre- 
quencies, this  result  is  true  no  matter  what  the  relation  is  between 

the  cp  ' s ; hence  the  third  field  could  be  anything).  Since  the  time 
P 

averaged  value  of  the  square  of  the  incident  field  is  simply  I E"/2  , 
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the  efficiency,  in  the  plane  wave  approximation,  is  given  by: 
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If  we  define  the  normalized  auto-convolution  of  the  power  spectrum  of 


the  incident  field: 


n 


. ■ 


where  I is  the  average  intensity  of  the  incident  field,  this  can 
be  rewritten  as: 


This  expression  is  very  general  in  that  it  gives  the  lineshape  of  the 
two-photon-resonance  for  any  power  spectrum  with  independent  modes, 
both  in  the  Doppler  and  the  pressure  broadened  regime.  In  an  arbi- 
trary case,  (2.25)  can  be  evaluated  numerically  (see  Appendix  A . Two 
special  cases  are  of  particular  interest,  because  they  occur  often  in 
practice,  and  because  they  also  suggest  a new  way  of  measuring  line- 
widths  of  non-allowed  transitions.  They  are  the  cases  when  either 
pressure  or  Doppler  broadening  dominates,  while  the  mode  spacing  is 
much  smaller  than  and  the  laser  bandwidth  is  much  wider  than  the  atomic 
linewidth  of  the  case  being  considered.  We  will  now  derive  closed  form 
expressions  for  the  efficiency  in  both  these  cases,  assuming  a Gaussian 
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power  spectrum  for  the  incident  radiation: 
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where  B = (£n  2)  2 times  the  F.W.H.M.  of  the  laser  spectrum.  Using 
this  definition  one  easily  finds,  by  substitution  in  . : 
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a . Pressure  broadened  regime,  \ B 


On  resonance,  one  finds  in  this  limiting  case: 
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Substituting  . and  . in  the  summation  over  k appearing  in 
(2.25),  we  find  for  this  sum: 
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Since  the  laser  bandwidth  is  assumed  to  be  much  wider  than  the  atomic 
linewidth,  the  exponential  in  the  sum  can  be  set  equal  to  one.  Because 
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the  mode  spacing  is  much  less  than  the  linewidth,  the  sum  can  be  ap- 
proximated by  an  integral.  The  final  result  is: 
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7T 

which  shows  an  inverse  linear  dependence  on  both  the  laser  bandwidth 
and  the  pressure  broadened  atomic  linewidth. 

b.  Doppler  broadened  regime,  A 7 ^ ' b 
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In  this  limiting  case,  on  resonance: 
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Again  we  have  to  evaluate  the  sum  over  k appearing  in  . . S . Using 
(2.27)  and  . . ' 1 ' : 
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Because  the  laser  is  assumed  to  be  much  wider  than  the  Doppler  width, 
the  term  involving  the  laser  width  can  be  neglected  in  the  exponential. 


Because  the  mode  spacing  is  much  less  than  the  Doppler  width,  the  sum 


can  be  approximated  by  an  integral.  The  result  is: 
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For  the  same  reasons  as  above,  this  can  be  approximated  by: 
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The  value  of  the  integral  in  this  expression  is  .5.  So  the  final  re- 


sult is  : 


I — 1 


Substituting  .33  and  . in  . we  find  for  the  efficiency 


- UO  - 


a*,., 


in  this  case: 


S = 


■ - A yj~  f - yjT  ) 


o !i  p p p 

(>V^Lcrav 

V 


(Z^kL/2)  (2.37) 


which  shows,  similar  to  the  pressure  broadened  case,  an  inverse  lin- 
ear dependence  on  the  laser  bandwidth  and  the  Doppler  broadened  atomic 
linewidth . 


Remarks 

The  inverse  linear  dependence  of  the  efficiency  on  $ and  7 
(2.30),  and  S and  7^  (2.37)  can  easily  be  arrived  at  using  simple 

arguments.  Assume  that  the  incident  field  consists  of  N modes  of 
equal  amplitude  E . In  the  pressure  broadened  regime,  the  response 
at  2 cUq  , and  hence  the  electric  field  generated  by  a pair  of  modes 
(m, n)  is  proportional  to  E ‘ /y  . Since  the  contributions  of  dif- 
ferent pairs  have  independent  phases,  powers  have  to  be  added.  So  we 
just  have  to  multiply  the  power  for  one  pair  of  modes  by  the  number  of 
pairs.  This  number  can  easily  be  found  as  follows.  If  we  keep  the 
mode  used  in  the  first  perturbation  step  fixed,  then  there  are  7 /Vi 
(for  An  < 7 < ?>  ) possibilities  for  the  mode  in  the  second  step.  This 
can  be  done  using  any  of  the  N modes  in  the  first  step.  Thus  the 
total  power  and  hence  the  efficiency  is  proportional  to: 


, this  can  be  written  as: 


; 

In  terms  of  the  average  intensity 

e ec 

In  the  Doppler  broadened  regime,  the  atoms  can  be  divided  into  packets 
of  width  7 . All  the  atoms  in  one  packet  are  mainly  driven  by  the 

same  set  of  modes,  while  atoms  in  different  packets  are  driven  by  in- 
dependent sets  of  modes  . Hence  in  this  case  the  total  radiated  power 
will  be  proportional  to  the  number  of  packets,  while  in  the  pressure 
broadened  regime,  where  all  atoms  respond  identically,  it  varies  as 
the  square  of  the  total  number  density.  This  means  that  compared  to 
the  other  case,  the  efficiency  is  down  by  the  number  of  packets,  i.e.: 

I" 

g «c  -22 
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Note  that  in  this  last  derivation  the  result  of  the  pressure  broadened 
regime  was  used,  which  means  that  it  is  only  valid  if  An  • • , not 
just  An  < 7^  . In  the  exact  derivation  we  did  not  have  to  make  this 
stronger  assumption.  It  should  be  noted  however  that,  if  the  mode 
spacing  is  larger  than  the  linewidth  of  individual  atoms,  we  have  a 
non-homogeneous  situation  where  certain  atoms  are  driven  much  harder 
than  others.  The  result  (2.37)  gives  the  efficiency  for  the  total 
ensemble  of  atoms.  In  using  it  under  non-homogeneous  conditions  one 
should  check  however  whether  the  response  of  strongly  driven  atoms  is 


I « NE2 
av 


av 
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not  saturating  (a  similar  remark  can  be  made  for  two-photon  absorption, 
where  bleaching  can  occur  . 

If  the  laser  spectrum  consists  of  only  a few  modes,  then  the  gen- 
eral expression  for  the  efficiency  (2.25)  should  be  looked  at  as  an 
ensemble  average,  i.e.,  it  represents  the  average  efficiency  for  a 
large  number  of  laser  shots. 

. . Time -Bandwidth -Limited  Pulses 

We  calculate  the  tripling  efficiency  for  an  electric  field  of  the 

form: 


E(t) 


-t 


COS 


V 


(2.38) 


Substituting  this  expression  in  (2.1): 


E (CD) 


EqT  I - [(a)-  no)'  t /U] 

1 e 

u ( 


-[(to  * 2T5  A I I 


The  polarization  set  up  at  the  third  harmonic  frequency  is  found  by 


substituting  (2.39)  into  (2.2): 


P(3)(t)  = _hK1E3(t)  X 


(JoL) 

U ^ > 
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do),  dtoo  e 


-(t2A)((<»1-'«0)2  + (®2-'»0)2l  -!(«>!  +°>2)t 
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_(T2/10[(a>1+a>o)2  + K,+a>Q)2]  -i(«1+»2)t  j 

/V  d°ldn2  e " ! 

JJ  (x>2Q  + + a>2  + i.7/2  j 


The  two  integrals  can  be  greatly  simplified  by  making  two  successive 
changes  of  variables: 


First  Change 


First  Integral 


CO, 

- CO 

1 

0 

m 

- . 

2 

0 

Second 


x = 

y = 


Integral 


«L  + “O 


CD  + 05 


Second  Change 
u = x + y 
v = x - y 
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The  result  is: 
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(2.-1 


The  integral  with  respect  to  v is  trivial,  while  the  ones  with  re- 
spect to  u can  be  quite  easily  done  if  looked  at  as  Fourier  trans- 
it 

forms  of  products.  The  result  is: 


p1' ' 1 (ti 


- 2KiEoE3(t) 


X ' le 


I 1 


- t : 


-?t/2  -i  ' t 

H(t)  e 


+ c .c , 


(2.142) 


where  H(t)  is  the  unit  step function  and  the  star  denotes  convolution. 
One  can  show  that  in  the  rotating  wave  approximation,  this  expression 


Since  we  are  dealing  with  a pulse  with  a certain  profile  in  time,  the 
most  meaningful  definition  for  the  efficiency  in  this  case  is  as  the 
ratio  of  the  energy  in  the  third  harmonic  pulse  to  the  energy  in  the 
fundamental  pulse.  In  the  plane  wave  approximation,  where  all  beams 
have  the  same  area: 


. V 

/ [E  ^ (t) r dt 
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•00 


00 
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dt 
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32  y[? 

T?  ^rT 


vn)2v?Lc 
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(j/A)^  X sin"  (2kL/2) 


where 


-00 


(2.Ut3 

The  quantity  is  a dimensionless  quantity,  while  J A is  the  en- 

ergy density  in  the  incident  pulse.  The  expression  (2.^tb)  gives  the 
shape  of  the  resonance  profile,  assuming  no  Doppler  broadening  (see 
later).  On  resonance,  in  the  usual  case  that  the  pulselength  is  much 
than  the  atomic  dephasing  time,  the  function  C ( t reduces  to 
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shorter 


the  probability  integral,  and  D^  can  easily  be  calculated  numerically 
giving  the  result  = O.o7  . Hence,  on  resonance: 

? = -55  x (J/A)r  x sin'  (AkL/2)  (2.h7 

As  mentioned  before,  this  result  neglects  Doppler  broadening.  For  the 
case  where  the  laser  bandwidth  is  considerably  wider  than  the  Doppler 
width,  which  is  usually  the  case  for  pulses  generated  by  a mode-locked 
laser,  it  can  simply  be  argued  however  that  Doppler  broadening  is  not 
going  to  affect  the  result  (2.^7).  The  reason  is  that  in  the  mode- 
locked  case,  in  contrast  to  the  case  of  independent  modes,  all  the 
phases  are  zero,  and  hence  all  the  atoms  are  excited  in  phase,  even 
though  different  packets  are  driven  by  different  bands  of  frequencies. 

Remarks 

The  dependence  of  (2.^6)  on  laser  linewidth  can  easily  be  under- 
stood as  follows.  Applying  a short  optical  pulse  with  center  frequency 
'ti  is  by  definition  equivalent  to  applying  a monochromatic  wave  of 
frequency  f'V  for  a time  equal  to  the  pulselength  T . The  polar- 
ization set  up  by  a monochromatic  wave  on  resonance  is  proportional  to 
E"* fy  . However,  the  build-up  time  on  resonance  is  T,  , and  since 
we  assume  that  this  time  is  much  longer  than  the  pulselength,  the 
polarization  will  only  build  up  to  a fraction  t/t  of  the  value 
for  an  infinite  wave.  Thus  the  value  of  the  polarization  for  a short 
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pulse  is: 


*5 


P =c 


X « E^  T 

T2 


The  energy  in  the  polarization  pulse  and  hence  in  the  radiated  elec- 


3 2 

trie  field  is  then  proportional  to  (E  f)  X T . Hence  the  efficiency 


is  proportional  to: 


EoT? 

F CC  OC  (j/A) 
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D.  TWO -PHOTON  ABSORPTION 


1 . Electric  Field  With  a Discrete  Spectrum  of  Independent  Modes 


We  calculate  the  two-photon  absorption  probability  per  second  for 
an  electric  field  of  the  form  (2.13),  using  (2.8).  First  we  calculate 
the  back -polarization  at  the  fundamental  frequency,  using  (2.3).  The 
result  is  an  expression  very  similar  to  (2.18),  except  that  the  tem- 


poral variation  is  at  o>  , and  the  constant  K1  gets  replaced  by  a 


1 


constant  K0  , where: 


N 


X 


^01^12 
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Hence : 
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Using  this  expression  we  van  calculate  PE  . Keeping  only  the  slowly 
varying  terms,  which  are  the  only  ones  that  will  contribute  to  the  time 
averaged  absorption  rate,  we  find: 


K. 


PE  = 


{ 


[-  FR(k)  sin  \J/  + F^(k)  cos  ip] 


2 V 27T7d  ni,  n,p,q 


where 


t = (m  + n - p - q)Aot  - (cp  + cp  - cp  - cp  ) (2.-9) 

Y ' r n m n p q 


According  to  (2.8),  we‘  now  have  to  calculate  the  time  averaged  value  of 
this  quantity,  i.e.: 


PE) 


lim 
T-»  » 


k f » 


dt 


We  find  immediately  the  condition 


'm  + n = p + q)  , and  hence: 


(PE) 
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m, n, p, q 
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E E E E FT  k cos  + cp 
m n p q 1 m n 
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Because  of  the  independence  of  the  phases,  the  only  terms  which  do  not 

average  out  in  this  quadruple  sum  are  those  for  which  -t  ' - c - 

m n p q 

such  that  the  cosine  factor  is  unity.  This  is  only  possible  if: 


( - ’p  l ‘ 

) °r  I 

| 9 = Cp  ( CT  c. 

\ n q V n p 


which  says  that  the  same  pair  of  modes  has  to  be  used  in  going  from  dc 
to  frequency  2c and  in  coming  back  down  from  205  to  dc . The  ex- 
pression : ,5!'  then  simplifies  to  a double  sum: 


(PE) 


K 
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E E;Ek-„Fi(k' 


(2.51 


This  expression  can  be  rewritten  in  terms  of  the  auto-convolution  func- 
tion G(k'i  , and  substituted  into  , . . Thij  gives  the  final  general 
resul t : 


2K^q  I2 

1 

(2)  2 0 av 

Wv 1 = X - 
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- SI 


As  it  should  be,  this  expression  only  contains  the  coefficients  k) 

of  the  quadrature  component  of  the  back -polarization . We  will  now 
again,  as  in  the  previous  section,  derive  closed  form  expressions  for 
the  special  cases  where  the  laser  bandwidth  is  much  larger  than  and 
the  mode  spacing  is  much  smaller  than  the  atomic  linewidth,  and  either 
pressure  or  Doppler  broadening  dominates. 


a . Pressure  broadened  regime,  .V  ■ y • g 
On  resonance: 


?/2 


(kAo)2  + (7/2) 2 


This  expression,  together  with  (2.27),  can  be  used  to  evaluate  the  sum- 
mation over  k appearing  in  (2.57  . Using  identical  arguments  as  in 
Section  C,  one  finds: 


/ , G(k)FI(k)  - . 

^2n72D  k * 

Substituting  this  expression  and  the  expression  for  K , we  find  for 

c 

the  two-photon  absorption  rate: 


W 


(2) 
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Note  that  this  expression  does  not  contain  the  atomic  linewidth. 


U 


b.  Doppler  broadened  regime,  .'V  • > • 5 


On  resonance: 


Fj (k)  = 


- k Sc  ' ' ■ 


Combining  this  expression  with  (2.27),  and  using  identical  arguments 
as  in  Section  C,  one  finds: 


y c(k)FI(k) 


1 2717^  k 


This  result  is  identical  to  that  of  the  pressure  broadened  case  (2.5k  , 

(2) 

and  thus  Wv“  is  given  by  (2.55  • 


Remarks 


As  for  the  efficiency  expressions  in  Section  C,  the  inverse  lin- 
ear dependence  of  the  two-photon  absorption  rate  on  laser  linewidth  can 
be  found  using  simple  arguments.  Assume  again  that  the  incident  field 
consists  of  N independent  modes  of  equal  amplitude  E.  The  time 
averaged  absorption  rate  (PE)  for  one  pair  of  modes  (m, n is  pro- 
portional to  E />  X E E */?  . Because  of  the  independence  of  the 

phases,  the  total  absorption  is  simply  found  by  multiplying  the  absorp- 
tion due  to  one  pair  of  modes  by  the  number  of  pairs  such  that  their 
sum  frequency  is  under  the  atomic  linewidth.  We  derived  in  Section  C 
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that  this  number  is  '//An  < N if  T ? ■ 


Hence : 


t '/ 

W ‘ « — — N 

7 An 


or,  in  terms  of  1 « NE'  : 

av 


W<2>  « ^ 


independent  of  the  atomic  linewidth.  In  the  Doppler  broadened  r«  • i:  t . 
since  absorption  is  an  "individual  atom"  phenomenon,  the  abou  expres- 
sion is  clearly  unchanged  if  the  laser  is  much  wider  than  the  Uppler 
width,  such  that  all  atoms  are  driven  equally  hard.  If,  in  the  Doppler 
regime,  the  mode  spacing  is  larger  than  the  linewidth  of  individual 
atoms,  a non-homogeneous  situation  arises.  One  should  cluck  in  that 
case  whether  the  atoms  which  are  driven  hardest  are  not  being  bleached. 

If  the  laser  spectrum  consists  of  only  a few  modes,  . must  be 
looked  at  as  an  ensemble  average. 


; . Time-Bandwidth-Limited  Pulses 


We  assume  an  electric  field  of  the  form  . , The  back-polar- 
ization is  then  given  by  an  expression  similar  to  . , but  with 

replaced  by  K , and  i replaced  by  : 


O r •. 

P ( t ) - K E e~l  : [C.,  t sin  t C t 

, l ' K l 


cos  t 


After  differentiation  and  multiplication  by  the  incident  field, 
find,  keeping  only  slowly  varying  quantities: 


we 


PE  - Kon  E e 


' t < T ^ 


Tht  total  probability  for  absorption  during  the  pulse  is  found  by 
substituting  this  expression  in  (2.9).  The  result  is: 


i D K J A 

W dt  = — 


where 


/ dv  e ' Cr(t\ 


The  quantity  D,  is  a dimensionless  quantity,  while  J /A  is  the  en- 
ergy density  in  the  pulse.  On  resonance,  for  the  usual  case  where  the 
pulselength  is  much  shorter  than  the  atomic  dephasing  time,  one  finds 
numerically  1).  . . Substituting  this  value  into  , to- 
gether with  tlie  expression  for  K . L}  , we  find  for  tht  total 

absorption  probability  on  resonance: 


W dt 


• '•  i", 


Note  that  this  expression  is  independent  of  the  atomic  linewidth.  In 


the  derivation  Doppler  broadening  was  not  considered,  but  as  in  Sec- 
tion C it  is  clear  that  the  result  : will  also  be  correct  in  the 

Doppler  broadened  regime,  if  the  laser  bandwidth  is  much  wider  than 
the  Doppler  width  such  that  all  atoms  are  driven  equally  hard. 

Remarks : 

The  dependence  of  .1  on  J A can  easily  be  understood  as 
folic ws.  As  explained  in  Section  C,  the  back-polarization  is  propor- 
tional te  E r Hence 

(2 ) 

W <x  PE  <*  E . T 

and  the  integrated  absorption  probability  will  be  proportional  to: 

CO 

J W!'2'  dt  ec  (eJt)t  cc  (j/A)C 


E.  QUADRATIC  KERR  EFFECT 

1 . Electric  Field  With  a Discrete  Spectrum  of  Independent  Modes 

We  will  calculate  the  Kerr  contribution  to  the  index  of  refrac- 
tion using  .1  and  .11  . Consider  a mode  with  frequency  • ml 
The  third  order  back-polarization  at  this  frequency,  due  to  the  mode 
itself  and  all  the  other  modes  is  given  by  an  expression  similar  t> 


(2.58),  but  with  p = n , i.e.: 


P!^(t;a)o  + mAx>)  = 
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V 


X > E 

-1  Z-rf  n 

2717  ‘D  n 


X ! F (k)  cos  [ (■  o.  + n*Qu>)  t - q>  ] 
/ K m 


+ FT(k)  sin  [ (uj  + nv\u)t  - Cf  ] 
I m 


The  component  in  phase  with  the  incident  field  can  be  written  as 


F (m  + n ' 

P;  . = - X V"*  E‘  ■ E It 

in  phase  2 / ^ n 1 771  m 

n V 27T7q 


Combining  this  expression  with  (2.1  ' and  (2.11),  we  find: 


I 


If  we  substitute  this  expression  together  with  the  definition  of 
(2.26)  into  (2.6U),  we  find: 
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If  we  let: 


nAn  = 


p P 


and  assume  that  the  mode  spacing  is  much  less  than  the  atomic  linewidth 
such  that  the  sum  can  be  approximated  by  an  integral,  this  can  be  re- 
written as: 


n ' 1 Kerr 


where 


t(m,S')  = / 


, , -x2/2B'2 

(x  - nv-xc  1 e dx 


(x  -toScd)2  + (7/2  V 


If  the  laser  bandwidth  is  much  wider  than  the  atomic  linewidth,  we  can 


use  the  asymptotic  form  of  F^(m,  ) 


(- 

(m\n)2  + (7/2 Y 


Combination  of  (2. o'  and  2/  i)  yields,  in  terms  of  the  average  in- 
tensity of  the  incident  field: 
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Finally,  substituting  the  definition  of  K 2.58  : 


'n  " 1 ' Kerr 
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the  absolute  value  of  (n  - 1)  is  maximum  when  the  detuning  aVn  ; ± . /2 
At  this  detuning  the  last  factor  in  (2.70)  equals  1/7 


b . Doppler  broadened  regime 

The  limiting  form  of  F^m  + n)  , on  resonance,  is  given  by  . .31  . 
The  summation  over  n in  (2.1  can  only  be  done  numerically  in  this 
case,  and  we  will  not  consider  it  in  detail  here. 


2 . Time-Bandwidth -Limited  Pulses 

The  in-phase  component  of  the  hack-polarization  is  given  by  . .5  : 


P t 


K K‘  C ’ t ' K t 


. 1 


39  * 


To  find  the  effect  of  this  polarization  on  the  propagation  of  the  pulse 
through  the  nonlinear  medium,  one  has  to  solve  the  T.W.E.  (2.12).  Note 
that  right  on  resonance,  the  function  C^(t>  vanishes,  and  the  effect 


ls  zero. 


Although  it  is  not  rigorous  for  short  pulses,  one  can,  analogous 
to  the  infinite  plane  wave  case,  calculate  an  index  of  refraction  from 
(2.  ’l),  using  (2.10)  and  (2.11).  The  result  is: 


n - 1) 

Kerr 


UnQK2  (j/a)  Cj.  (t) 


From  the  maximum  value  that  this  expression  obtains  during  the  pulse, 
one  can  calculate  a maximum  value  for  1c  X L , which  should  give  a 
good  indication  of  the  importance  of  the  effect. 


F.  MEASUREMENT  OF  THE  LINEWIDTH  OF  ELECTRIC  DIPOLE  FORBIDDEN 
TRANSITIONS 


Several  methods  are  currently  available  to  measure  the  linewidth 
of  non-allowed  transitions,  ^ but  the  required  laser  sources  are  not 
often  available.  Using  the  results  of  Section  C,  it  is  possible  to 
quickly  and  conveniently  measure  such  linewidths,  even  when  they  are 
considerably  narrower  than  the  Doppler  width.  The  method  requires  a 
laser  source  which  can  be  several  orders  of  magnitude  wider  than  the 


atomic  linewidth  being  measured. 
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We  will  first  discuss  the  theory  on  which  the  method  is  based,  and 

2 2 

then  give  experimental  results  for  the  broadening  of  the  Js  S - 3d  D 
transition  in  Na. 


1 . Theory  of  the  Measurement 


Because  of  the  dependence  of  the  efficiency  in  the  pressure  broad- 
ened regime  on  the  linewidth  y (2.30),  one  could  conceivably  measure 
this  linewidth  by  measuring  the  efficiency  in  one  coherence  length 
versus  the  pressure  of  the  nonlinear  species  (self-broadening),  or  the 
pressure  of  a foreign  gas  (foreign  gas  broadening).  This  would  require 
however  detailed  knowledge  of  atomic  parameters  (matrix  elements),  and 
also,  as  the  pressure  is  changed,  the  efficiency  may  be  influenced  by 
effects  like  molecular  absorption,  change  of  the  coherence  length 
by  the  foreign  gas,  etc.  The  requirement  to  know  the  atomic  parameters 
may  be  avoided,  and  those  interfering  effects  normalized  out  by  making 
efficiency  measurements  on-resonance  and  of f-resonance,  and  taking  a 
ratio.  Of f-resonance,  the  efficiency  is  independent  of  both  the  laser 
bandwidth  and  the  atomic  linewidth,  but  it  will  still  be  influenced  by 
the  other  processes  mentioned  above.  For  sufficient  detunings  (several 
laser  bandwidths),  the  efficiency  is  given  by  (2.25): 
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The  ratio  of  the  efficiency  on-resonance  to  the  efficiency  of f -resonance 
is  found  by  combining  (2 .30)  and  (2.73): 


R = 2 


(cr>  - 2o)  V 
v 20  0; 


(2.74) 


7T  + 4 


. 


One  can  thus  determine  the  pressure  broadened  linewidth  y from  the 
measured  \ alues  of  R and  R : 


y = 1.12  X y X — 

R 


The  detuning  and  the  laser  bandwidth  in  this  expression  can  further  be 
eliminated  by  normalizing  to  the  low  pressure  limit  of  R which  is 
found  by  combining  (2.37)  and  (2.73)-  The  result  is: 


where  is  given  by  .T  . This  procedure  does  not  require  de- 

tailed knowledge  of  the  laser  spectral  composition,  but  is  only  valid 
in  the  pressure  broadened  regime.  Practically,  this  is  usually  the 
only  region  where  values  of  7 are  of  interest,  and  extrapolation  to 
lower  pressures  can  usually  be  made  with  confidence. 

A more  general  determination  of  y , valid  at  all  pressures,  can 
be  made  using  an  exact  numerical  evaluation  of  the  expression  for  the 


1 

efficiency  (2.25).  Using  the  fixed  parameters  7 , & , and  Zto  , the 

ratio  R /R  - = £ /£  ^ is  computed  as  a function  of  7 and  plotted 
c cO  c cO  r 

(the  subscript  c refers  to  calculated  values,  m to  measured 

values).  The  experimental  procedure  is  identical  to  that  explained 

above.  Two  sets  of  efficiency  measurements  (in  one  coherence  length', 

on-  and  off-resonance,  are  made  and  the  normalized  ratio  R /R  is 

m mO 

plotted  versus  the  pressure  of  the  broadening  agent.  The  normalized 
curves  of  measured  and  calculated  efficiency  are  then  used  to  deter- 
mine 7 from  the  relation  R (P)/R  = R (7)/r  n . Thus  for  any 

pressure  P , the  linewidth  7 may  be  read  from  the  point  on  the 

calculated  curve  associated  with  the  value  RPR 

m " mO 

The  two  procedures  are  illustrated  below. 

2 2 

2 . Measurement  of  the  Broadening  of  the  Na  3s  S - 3d  D Transition 

2 2 

The  foreign  gas  broadening  by  He  of  the  Js  S - Jd  D transition  in 
Na  was  measured  using  the  above  techniques. 

A pumping  wavelength  of  6856  % was  obtained  from  an  optical  para- 
metric oscillator,  and  the  idler  frequency  at  2.3"  p was  used  as  the 

r 

input  frequency  05  to  generate  the  observed  lower  sideband  signal  at 
2cOq  - 05  (U006  X)  . The  pumping  radiation  had  a bandwidth  of  2 cm  ' 
and  a mode  spacing  of  0.028  cm  ^ . Both  incident  beams  had  a peak  power 
of  about  1 kW  and  were  focused  into  a small  Na  heat  pipe  at  a tempera- 
ture of  about  700°K,  yielding  a Doppler  linewidth  7^  - . cm"'  f 

and  a Na  pressure  of  0.75  torr  (1  x 10^  cm 


Efficiency  measurements  were  made  both  on-  and  of f-resonance 


while  varying  the  He  pressure  from  a few  torr  to  two  atmospheres . 
Figure  1 shows  the  ratio  R^P)/!!  as  a function  of  the  He  pres- 
sure, and  Fig.  2 is  a plot  of  the  computed  ratio  Rc(7)/R,  for 
this  case.  Figure  1 also  shows  the  atomic  linewidth  y as  a func- 
tion of  pressure,  both  using  the  exact  method  and,  for  P torr  , 

using  the  approximate  method.  Note  that  the  two  techniques  are  con- 
sistent, and  that  the  numerical  technique  permits  one  to  measure 
pressure  broadened  linewidths  several  times  smaller  than  the  Doppler 
width.  The  measured  broadening  coefficient  is  .J  cm  * per  atmosphere 
of  He . 

A similar  procedure  was  used  to  determine  the  sel f-broadening  of 
this  transition  by  varying  the  Na  density  from  . 5 torr  to  150  torr. 
Unfortunately,  at  these  high  Na  densities,  it  was  necessary  to  intro- 
duce 300  torr  of  He  buffer  gas  to  protect  the  cell  windows,  and  the  Na 
self-broadening  was  not  observable  above  the  0.1  cm  ^ He  broadening. 
Thus  we  infer  a self-broadening  of  less  than  0.1  cm  ' at  a density  of 
1018  cm'5. 

From  a device  point  of  view  these  results  are  very  important. 

They  indicate  that  a Na  up-converter  (operating  aroupd  ■ may  be 
operated  at  pressures  of  about  100  torr  without  a significant  reduc- 
tion in  efficiency  due  to  broadening  of  the  atomic  linewidth,  and 
second,  that  it  is  important  to  keep  the  buffer  gas  pressure  as  low 
as  possible.  Figure  1 shows,  for  example,  that  at  a He  buffer  gas 
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Fig.  2--Calculated  ratio  R 'R  versus  atomic  linewidth. 


pressure  of  about  300  torr,  the  efficiency  is  about  a factor  of  three 
less  than  the  value  with  no  He  at  all. 
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CHAPTER  III 

INFRARED  IMAGE  UP-CONVERSION  IN  ALKALI  METAL  VAPORS 
A.  INTRODUCTION 

The  interest  in  frequency  shifting  radiation  from  the  IR  into  the 
visible  or  near  UV  stems  mainly  from  two  factors.  First,  there  are  now 
available  very  efficient  low-noise  detectors  of  visible  radiation  such 
as  the  human  eye  and  photomultipliers,  whereas  most  infrared  detectors 
are  inefficient  (and  expensive)  by  comparison.  A second  attractive 
feature  of  this  method  of  detecting  infrared  is  that  the  up-converter 
and  the  visible  detector  following  it  can  all  operate  at  room  tempera- 
ture, whereas  the  competitive  detectors  in  the  infrared  commonly  oper- 
ate between  ^.2°K  and  77°K. 

In  this  chapter,  infrared  imaging  using  two-photon-resonant  fre- 
quency converters  is  discussed  in  detail.  Both  passive  and  active 
imaging  are  considered,  but  the  emphasis  is  on  active  imaging.  The 
reason  for  this  is  that  a recent  study*  has  shown  that  active  imaging 

■v 

systems  using  such  an  up-converter  look  attractive  when  compared  to 
competing  techniques  using  photoconductive  detectors  or  quantum  coun- 
ters, but  they  compare  rather  poorly  with  these  other  techniques  for 
passive  imaging  of  thermal  radiation.  The  main  reason  behind  this  is 
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that  the  up-converters  have  to  use  a pulsed  laser  for  a pump,  resulting 
in  a low  duty  factor.  In  the  case  of  active  imaging,  the  IR  radiation 
itself  is  pulsed,  and  this  duty  factor  problem  is  much  less  severe. 

We  will  start  with  a brief  discussion  of  the  performance  criteria 
for  passive  and  active  imaging  systems.  This  will  be  followed  by  a 
description  of  the  various  optical  systems  that  can  be  used  with  up- 
converters  . The  next  section  deals  with  field  of  view,  resolution 
and  bandwidth.  Then  follows  the  main  section  of  this  chapter,  which 
is  a detailed  study  of  the  efficiency  of  such  up-converters,  the 
physical  processes  limiting  it,  optimization  procedures  and  a com- 
parison with  crystal  up-converters.  The  last  section  is  a discussion 
of  the  experimental  results.  It  includes  a description  of  the  system 
that  was  used  to  obtain  the  first  images,  together  with  a comparison 
between  measured  and  calculated  efficiencies. 

Whenever  possible,  theoretical  results  will  be  applied  to  the 
example  of  2.*ph  p up-conversion  in  cesium  vapor,  using  the  process 
1.0790  p + 1.0790  p + 2.0^  p H . The  interest  in  this  par- 

ticular system  stems  from  the  fact  that  the  high  power  Nd:BeL  solid 
state  laser  at  1.0790  p has  an  accidental  two-photon-coincidence  with 
the  6s-7s  transition  in  cesium. 

B.  INFRARED  IMAGING  SYSTEMS 

Performance  Criteria  for  Iinapiiii;  Systems 


1. 


In  passive  imaging,  an  object  is  detected  by  the  thermal  radiation 
it  emits  at  any  temperature  above  absolute  zero.  Different  parts  of 
the  object  which  are  at  slightly  different  temperatures  emit  different 
amounts  of  radiation.  The  sensitivity  of  a detector  or  an  imaging  de- 
vice is  evidently  related  to  the  relative  magnitudes  of  its  response 
to  a change  AT  , and  the  random  fluctuations  in  its  response  due  to 
various  noise  processes.  Therefore,  one  'widely  used'  performance  in- 
dicator of  a thermal  detector  is  the  noise-equivalent  temperature  dif- 
i ference  (NETD),  which  may  be  defined  as  that  effective  temperature 

I difference  that  causes  a response  equal  to  the  rms  noise  fluctuations 

in  the  response  of  the  detector,  i e.,  giving  a signal-to-noise  ratio 
of  unity.  We  will  now  derive  an  expression  for  the  NETD. 

Hie  number  of  photons  received  per  second  by  a detector  with  area 
A and  field  of  view  All  , looking  at  an  object  of  emissivity  <:  (> 
and  temperature  T , is  given  by:^ 


In  this  expression  p is  a factor  that  takes  into  account  the  response 
of  the  detector  for  orthogonal  polarizations.  If  botli  polarizations  are 
up-converted  with  the  same  efficiency,  p equals  1 . If  only  one  polar- 
ization is  up-converted,  its  value  is  1 . If  both  polarizations  arc  up- 
converted  with  different  efficiencies,  which  is  the  case  for  two-photon- 
resonant  frequency  converters  making  use  of  a s to  d transition,  one 
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has  1 < p < 2 . The  dependence  of  third-order  processes  on  polar- 

ization is  discussed  in  detail  in  Appendix  E.  Since  the  bandwidth 
of  vip -converters  is  narrow  compared  to  the  center  frequency  of  the 
up-converted  band,  the  integral  in  (3.1'  can  be  approximated  by  the 
product  of  the  integrand,  evaluated  at  this  center  frequency,  and  the 
bandwidth . Thus 


pt  (X^  ) Af  AAd 


hc/X ,kT 

l 

e 

where  X^  is  the  wavelength  at  the 
the  bandwidth  of  the  up-converter, 
tion  element  is  found  by  setting  the 
equal  to  the  diffraction  angle,  i.e. 
of  photons  per  resolution  element  is 


- 1 XT 

l 

center  of  the  band,  and  Xf  is 
The  number  of  photons  per  resolu- 
solid  angle  in  this  expression 

p 

, Aid  XT  A . Hence  the  number 

i 

approximately  given  by 


n 


pe(Xi)  Af 

hc/x . kT 
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where  the  temperature  T is  the  average  temperature  of  that  element. 
Starting  from  this  expression,  we  can  derive  an  expression  for  the 
NETD.  The  derivation  presented  here  follows  largely  Falk  and  Tiffany.^ 
Adjacent  resolution  elements  of  the  object  have  slightly  different  tem- 
peratures, and  the  "signal"  is  given  by  the  difference  between  the  num- 
ber of  photons  received  from  two  adjacent  spots  with  a temperature  dif- 
ference \T  . Consider  two  such  spots,  emitting  at  average  rates  n^ 
and  n . The  signal  received  at  the  input  of  the  up-converter,  during 


h. 


"l 


a time  t 


is  then  equal  to: 


S. 

in 


(3-M 


For  Poisson  noise,  due  to  fluctuations  in  time  of  arrival,  the  noise 
is  equal  to  the  square  root  of  the  average  number  of  received  photons: 

Nin  = t(nx  + njt];  (3. 5'1 

From  the  expressions  for  the  signal  and  noise  at  the  input,  the  signal 
and  noise  at  the  output  can  be  calculated.  We  assume  up-conversion 
with  piioton  efficiency  , followed  by  detection  using  a device 

with  efficiency  q , both  processes  being  noisefree.  Thus: 
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We  will  use  this  expression  now  to  calculate  the  NETD,  where  the  dif- 
ference np  - n,  is  due  to  a tempera,  .re  difference.  One  could  also 
use  this  expression,  however,  to  calculate  the  S N ratio  for  detection 


of  a difference  in  reflectivity  between  adjacent  spots  which  could  be 


used  as  a performance  criterion  for  active  imaging  . Using  for  nj 
and  n.  expression  (j.j),  evaluated  at  temperature  T and  T • -T  . 
we  find: 
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By  definition,  the  NETD  is  a AT 
noise  ratio  at  the  output  of  the 


NETD  = 


(j_ 

pe(\)  Af  ' 

i>r 
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such  that  the  resulting  signal-to- 
imaging  system  equals  unity.  Hence: 
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Example : 

We  will  use  typical  parameters  for  a Na  up-converter,  operating 
around  9 Mn  the  p -»  12  p atmospheric  window. 

e(X.)  = 1 

T = JOO°K 

f.  = 1100  cm'1 

1 

Af  = 200  cm'1 

ePH  = °*25 

^ - 0.3 

p - 1.25 

t 1 psec 

— > NETD  1.5°K 


b. 


Active  imaging 


For  the  case  of  active  imaging,  where  the  object  is  illuminated 
using  a high  power  laser  and  detected  by  the  radiation  it  reflects  to- 
ward the  up-converter,  the  NETD  is  clearly  not  a useful  performance 
criterion.  One  possible  application  of  active  imaging  might  be  where 
one  is  trying  to  detect  an  object  with  some  finite  reflectivity  in  a 
non-reflecting  ("dark"  background.  A reasonable  performance  crite- 
rion for  this  case  is  the  signal-to-noise  ratio  for  the  detection  of 
the  difference  between  a spot  in  the  "black"  background,  and  a spot 
on  the  reflecting  object,  e.g.,  spots  1 and  2 in  Fig.  jj . A "good" 

S/N  ratio  would  then  permit  one  to  tell  the  shape  of  the  object.  We 
will  use  a constant  reflectivity  Ref  for  the  entire  object.  If 
this  is  not  the  case,  one  could  use  as  one  possible  choice  the  mini- 
mum value  of  Ref  , in  which  case  one  would  be  able  to  see  some  de- 
tail in  the  object. 

The  noise  sources  to  be  considered  in  this  case  are  the  shot 
noise  due  to  the  signal  itself,  thermal  radiation  coming  from  the 
object  because  it  has  a non-zero  temperature,  and  dark-current  noise 
in  the  phototube  following  the  up-converter.  We  will  derive  an  ex- 
pression for  the  signal-to-noise  ratio,  and  we  will^ also  find  a lower 
limit  for  the  power  of  the  illuminating  laser  required  in  order  to 
make  the  reflected  radiation  stronger  than  the  thermal  radiation. 


The  geometry  for  active  imaging  is  shown  in  Fig.  it.  We  will 


assume : 


SI  < Si  < 

o upc  ~ L 


If  we  call  T)  the  overall  system  efficiency,  we  have  (Fig.  3) 
Sout  = (ni  " "2)  X n X t = (ns  ~ 0)  X t|  x t 
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where 

i , = dark  current 

d 

t = integration  time 

ng  - number  of  reflected  photons  per  sec 
per  resolution  element  on  the  object 

n^  = number  of  thermal  photons  per  sec 

per  resolution  element  on  the  object 


(S/N 


out 


+ n. 


2(id/e)] 


X (h  X t) 


Wc  will  now  calculate  the  number  of  photons  reflected  by  the  object 
per  sec  per  resolution  element,  and  also  the  number  of  thermal  pho- 
tons per  sec  per  resolution  element.  If  the  illuminating  laser  emits 

N photons  per  second  into  a solid  angle  12  , then  the  number  of 

L •*-» 

reflected  photons  received  per  second  is: 

12  Ref  1 7TD' 

N = Nt  x — X X — r X (5.12) 

S L d ) 

°L  * r 

In  this  expression  D is  the  diameter  of  the  up-converter,  and  use 

was  made  of  the  fact  that  2 > 12  .If  the  total  number  of 

upc  obj 

resolution  elements  of  the  up-converter,  corresponding  to  the  full 

solid  angle  12  ^ > is  equal  to  R , then  the  number  of  resolution 

elements  on  the  object,  which  spans  a solid  angle  12  , equals 

R X 12  / 12  . Hence  the  number  of  reflected  photons  per  resolu- 

o upc 

tion  element  (on  the  object  is: 

N 12  /R  Ref  1 7tD 

n = 2 = N X — — X X — X (3.I3) 

S R x 12  /l2  12T  7T  r'  k 

o u L 

The  total  number  of  thermal  photons  emitted  by  the  object  at  tempera- 
ture T equals: 

p ttD2A 

Nb  = X Af  x (12  x r")  x (3.  A) 

r‘ 

where  J,  is  the  number  of  photons  emitted  per  s"c  per  unit  area, 
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is  given  by 


per  solid  angle,  per  unit  bandwidth.  Using  (3 .2] , 

Pt(\) 

JX  = _ hc/X.kT 

^(e  1 -1) 

This  radiation  comes  also  from  R X Q '.i  resolution  elements,  and 

o u 

thus  the  number  of  photons  per  resolution  element  is  equal  to: 


flu  TTD" 

n,  = J,  X Af  x — x 

b - „ , 

R It 


By  combining  (3.13  and  (3.16  , we  find  the  minimum  number  of  laser 
photons  per  second  required  to  make  the  reflected  radiation  dominate. 
The  result  is: 


I X X Af  x (U^r 


If  the  laser  power  satisfies  this  condition,  the  detection  is  called 
signal-shot-noise-limited  (neglecting  the  dark  current  . If  we  sub- 
stitute (5.13)  and  (3 .16)  in  the  expression  for  the  signal-to-noise 
ratio  (3.11  , and  neglect  the  dark  current  contribution,  we  find  in 
the  signal-shot-noise-limited  regime: 


(S/N)0ut  = (t>s  X tj  x t.)  = - D 
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Example  (Cs  up-converter  at  £.9^  p. 


12  = 1 X 10 

u 


D = 10  cm 


R = 10 


iL  = il 

L u 


Ref  = .01 


■n  = 0.1 


t = 10  nsec 


e(\) 


2 .9  1 P 


100  cm  (blue  filter  behind  up-converter 


°K 


P 


one  finds : 


yj' 


To  obtain  a signal -to-noise  ratio  of  unity,  the  required  laser  power  is 
: 'j50  kW.  Using  (3 . 1 for  the  reflected  radiation  to  dominate,  the 


required  power  is  about  200  times  less,  and  thus  we  are  indeed  in  the 
signal-shot-noise-limited  regime  . 


2 . Optical  Systems  for  Up-Converters 

We  will  only  briefly  discuss  the  various  optical  systems  that  can 

be  used  with  image  up-converters.  For  a more  complete  discussion,  see 

Refs.  16  and  18 . Three  basically  different  types  of  optical  systems 

18 

for  up-conversion  have  been  discussed  by  Andrews ; they  are  il lus trated 
in  Fig.  5. 

In  a Type  I optical  system,  the  object  is  imaged  in  the  up-conver- 
ter, and  the  up-converted  radiation  is  then  imaged  onto  a detector.  The 
image  of  a point  on  the  IR  object  is  formed  in  this  case  by  a spherical 
wavefront  converging  on  a point  in  the  up -converter . Up-conversion  em- 
ploying a Type  I optical  system  is  also  called  up-conversion  in  image 
space . 

In  a Type  II  optical  system,  nearly  plane  waves  from  the  IR  object 
are  up-converted,  and  the  up-converted  plane  waves  are  imaged  onto  the 
detector . 

In  a Type  III  optical  system,  the  object  and  the  up-converter  are  1 

placed  respectively  in  the  front  and  hack  focal  plane  of  a lens.  IR 
radiation  from  different  object  points  gives  plane  waves  with  different 
directions  in  the  up-converter.  The  up-converter  and  the  detector  are 

• ! 

respectively  in  the  front  and  back  focal  planes  of  a second  lens,  and 
this  lens  transforms  every  up-converted  plane  wave  into  a point  in  the 
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5__Xypes  of  optical  systems  for  IR  up-conversion. 


detector  plane.  Up-conversion  employing  a Type  III  optical  system  is 


also  called  up-conversion  in  Fourier  space. 

Because  with  a Type  II  system  there  is  no  freedom  in  matching  the 

incident  IR  radiation  to  the  area  and  field  of  view  of  the  up-converter, 

this  type  is  of  very  limited  use.  The  resolution  of  a Type  I system  is 

limited  below  the  diffraction  limit  by  the  length  of  the  medium,  while 

18 

this  is  not  the  case  for  a Type  III  system.  In  Appendix  B it  is  shown 
in  detail  that  the  resolution  of  a Type  III  system  is  only  limited  by 
diffraction,  assuming  monochromatic  light.  The  magnification  of  a 
Type  III  system  is  also  derived  there.  All  types  suffer  from  chromatic 
aberration.  For  a Type  III  system  this  is  discussed  in  Appendix  C. 
According  to  Ref.  It,  this  aberration  might  in  principle  be  corrected 
for  in  a Type  III  system,  which  would  probably  be  required  in  the  case 
of  broadband  thermal  up-conversion.  For  active  imaging  chromatic  aber- 
ration is  of  little  importance  since  the  radiation  is  very  narrovoband . 

In  a Type  I system,  the  output  is  the  product  of  the  IR  image  and  the 
square  of  the  pump  profile.  This  poses  very  severe  uniformity  problems 
on  the  pump  beam  for  the  case  of  low-contrast  imaging  thermal  imaging  . 
In  a Type  III  system,  on  the  other  hand,  each  point  of  the  object  gives 
rise  to  a plane  wave  which  is  up-converted  by  the  ent i re  pump  beam,  and 
as  a result  the  pump  beam  quality  is  much  less  important.  For  all  these 
reasons  a Type  III  system  is  usually  optimal.' 


C.  FIELD  OF  VIEW,  RESOLUTION,  AND  BANDWIDTH.  APPLICATION  TO  THE 

2.91-  n CESIUM  UP-CONVERTER 

1 . Active  Imaging 
a . Field  of  view 

Infrared  rays  incident  on  an  up-converter  at  different  angles  are 
up-converted  with  a different  k-vector  mismatch,  resulting  in  an  angle- 
dependent  efficiency.  Clearly  this  response  should  be  a "smooth"  func- 
tion, without  any  zeros.  There  are  two  possible  ways  for  obtaining  a 
ripple-free  response  versus  IR  angle.  A first  possibility  is  to  use  a 
cell  that  is  less  than  one  coherence  length  long.  At  the  high  pres- 
sures where  one  wants  to  operate  the  up-converters,  this  length  is  very 
short,  however,  typically  being  only  a small  fraction  of  a millimeter. 
This  is  usually  impractical,  both  because  of  the  mechanical  problems  oi 
building  such  a short  cell,  and  because  the  resulting  field  of  view  is 
so  large  that  it  becomes  hard  to  design  a cell  with  a large  enough 
"physical"  field  of  view.  The  second  possibility  is  to  make  the  medium 
reasonably  long  (there  are  also  limitations  there,  as  explained  in  Sec- 
tion D),  and  phasematch  over  its  entire  length.  We  will  discuss  only 
the  second  case. 

The  k-vector  mismatch  caused  by  the  angle  between  the  IK  rays  and 
the  direction  of  the  sum  of  the  two  pump  vectors  is  given  by: 

k.e2 
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The  efficiency  varies  versus  angle  according  to: 


(5.20) 


and  thus  the  solid  angle  of  the  up-converter  (k,_  points  is  given  by; 
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Note  that  this  solid  angle  is  exactly  equal  to  the  maximum  solid  angle 
as  allowed  by  the  difference  in  curvature  between  the  1R  and  visible 
wavefronts  in  the  up-converter,  as  discussed  in  Appendix  B.  By  using 
off-angle  phasematching,  ' where  a mismatch  .'ikL  - n is  allowed  for 
zero  1R  angle,  the  field  of  view  can  be  increased  by  a factor  of  two, 
but  in  view  of  the  results  of  Appendix  B this  may  be  unacceptable. 

b . Methods  of  phasematching 

Three  possible  techniques  for  phasematching  the  up-conversion 
process  will  be  discussed.  They  are  inert  gas  phasematching,  self- 
phasematching, and  phasematching  by  angle.  We  will  now  discuss  each 
of  these  techniques  in  detail. 

Phasematching  by  adding  an  inert  gas'  is  probably  not  a useful 
technique.  The  reason  is  that  because  of  the  high  pressure  ratios  re- 
quired, the  inert  gas  pressure  becomes  excessively  high  at  the  metal 
vapor  pressures  where  one  wants  to  operate  the  up-converter  on  t lie 


order  of  10  atmospheres).  This  leads  to  a large  broadening  of  the 
linewidth  of  the  two-photon  transition,  not  to  mention  the  practical 
problems . 

The  second  method  is  to  select  the  IR  wavelength  such  that  the 
medium  is  sel f-phasematched . This  method  is  the  simplest  one,  but  it 
requires  an  illuminating  laser  that  can  be  tuned  to  the  phasematching 
frequency.  We  will  now  derive  an  expression  for  the  detuning  from 
the  upper  level  where  the  process  is  phasematched,  by  setting  .'k 
equal  to  zero.  For  colinear  waves,  the  expression  for  the  k-vector 
mismatch  is: 
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The  factors  (n  - 1 in  this  expression  can  be  evaluated  using 

lo 

Sellmeier's  equation.  Assuming  no  population  in  the  s level,  the 
result  for  the  cesium  up-converter  is: 
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where  &n 
used  in  the 


is  in  wavenumbers,  and  only  the  6p  and  p doublets  were 
calculation.  If,  because  of  two-photon  absorption,  there 


is  some  population  in  the  two-photon-resonant  level,  it  will  have  only 
a small  effect  on  the  index  of  refraction  at  the  pump  and  the  sum 


frequencies,  but  it  will  considerably  affect  the  index  at  the  1R,  since 
the  IR  frequency  has  only  a small  detuning  for  the  transition  from  the 


resonant  level  to  the  upper  level  (7s-7p  in  cesium).  The  expression 
for  the  mismatch,  including  this  contribution,  is,  for  the  cesium  up- 
converter  : 
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where 


f = 100  X 


percentage  population 


N 


in  the  s level 


5. 


The  self-phasematching  frequency  is  found  by  setting  Ak  equal  to  zero. 
For  the  cesium  up-converter,  the  result  is: 


This  result  says  that,  as  the  population  of  the  7s  level  is  increased, 
the  phasematching  IR  frequency  will  originally  increase  such  that  the 
detuning  from  the  upper  level  becomes  smaller.  At  a smaller  detuning 
the  efficiency  is  higher,  but  absorption  botli  at  the  IR  and  the  sum 
frequency  has  to  be  checked  see  Section  D . As  the  fraction  f is 
increased  further,  the  process  will  eventually  phasematch  tor  negative 


detunings,  i.e.,  above  the  upper  level. 


The  third  technique  is  more  general  than  the  second  one  in  that 


Ik 


it  can  be  used  for  a range  of  IR  frequencies,  but  it  requires  a more 
complicated  system.  It  can  easily  be  understood  from  the  diagrams 
shown  in  Fig.  I.  If  the  sum  frequency  does  not  come  very  close  to 
the  upper  atomic  level,  the  dominant  contribution  to  the  k-vector 
mismatch  is  due  to  the  pump.  Since  the  frequency  of  the  pump  is  be- 
low the  first  doublet  of  the  metal  vapor,  it  sees  an  index  of  refrac- 
tion larger  than  1,  and  thus  the  vector  2k^  + k^  is  longer  than  the 
vector  kg  . By  splitting  the  pump  beam  in  two  beams  and  recombining 
them  under  a small,  angle,  the  mismatch  .^k  can  be  made  zero.  For 
small  angles  cp  , the  mismatch  caused  by  the  small  angle  . be- 
tween the  pump  beams  equals: 
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2<»  cp 

£k(cp)  - — 1 - X — (3  - t 
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By  adding  this  term  to  the  expression  the  total*  mismatch  for 

the  cesium  up-converter  becomes: 
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For  a given  IR  frequency  or,  equivalently,  a given  detuning  \ 
the  angle  cp  can  be  found  which  makes  this  expression  zero.  The 


(5,28) 

For  large  detunings  Ao_.  , the  angle  c.  approaches  an  asymptotic 
value.  In  Fig.  ' the  phasematching  angle  is  plotted  versus  detuning, 
for  a number  density  of  10^  cm  , and  for  various  fractions  f 

c . Resolution 

The  general  expression  for  the  resolution  of  an  up-converter  using 
a Type  III  optical  system  is  derived  in  Appendix  B.  If  we  substitute 
for  the  solid  angle  the  result  (5.21),  we  obtain: 

(3-29 


d. 


Relation  between  efficiency,  resolution,  and  number  densitv 


r 


Using  (3.29),  this  can  be  written  in  terms  of  the  number  of  resolution 
elements,  as: 


•- 


This  is  a key  relation  for  up-converters.  It  says  that,  for  a given 
efficiency  and  resolution,  the  required  power  is  inversely  proportional 
to  the  number  density  of  the  metal  vapor.  Thus  one  wants  to  operate  an 
up-converter  at  as  high  a pressure  as  possible  see  Section  D for  lim- 
iting effects ) . 


e.  Bandwidth 


The  efficiency  depends  on  the  detuning  from  the  upper  level  in  the 
following  way: 


(3  .31) 

where  Ak  is  given  by  (3.27  , and  also  depends  on  An,  . If  the  de- 
tuning is  changed  by  a small  amount  "(A:'..)  from  its  original  value 
An  , the  resulting  mismatch  is: 
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The  sine'  factor  in  the  efficiency  is  down  to  of  its  peak  value 

when  .'NkL  = 7T  . From  this  a maximum  value  for  5(Ai^)  can  be  cal- 
culated. We  will  illustrate  this  with  an  example  for  the  cesium  up- 
converter  . 


Example : 


We  find: 


2h  -5 

N = 10 m ' 


20  cm 


L = 2 X 10  " m 


f = 0 


= - O.65  x S(Ad  ) 

2 ^ 


and  thus 


[5(An  )] 

3 max 


~ 2 . 5 cm 


So  the  bandwidth  in  this  example  is  about  5 cm  • Note  that  the  fac- 
tor  l/(An_)  in  the  expression  for  the  efficiency  is  nearly  constant 
over  this  narrow  width.  In  general,  however,  it  also  affects  the  band- 


width, and  both  factors  should  be  taken  into  account . 


Passive  Imaging 


a.  Field  of  view  and  bandwidth 


For  passive  imaging  applications,  an  up-converter  should  have  a 
large  IR  bandwidth,  as  is  clear  from  the  expression  for  the  noise 
equivalent  temperature  difference  (3.9).  Because  of  the  large  band- 
width requirement,  phasematching  cannot  be  used  as  a means  of  getting 
a smooth  angular  response  in  this  case.  If  a cell  is  used  that  is 
several  coherence  lengths  long,  at  some  IR  frequency,  the  response 
vertus  angle  will  be  oscillatory  with  zeros  at  certain  angles.  How- 
ever, because  of  the  dependence  of  the  mismatch  on  the  detuning  \ , 

the  peaks  and  zeros  for  different  IR  frequencies  do  not  occur  at  the 
same  angles  As  a result,  a smooth  frequency  integrated  response 
versus  angle  can  still  be  obtained  by  using  a cell  that  is  several 
coherence  lengths  long  for  all  or  most  of  the  frequencies  in  the  band 
considered.  As  in  the  active  imaging  case,  one  has  the  possibility 
of  splitting  the  pump  beam  in  two  beams  at  a small  angle,  in  which 
case  one  has  a free  variable  that  can  be  chosen  for  an  optimum  re- 
sponse. The  problem  just  described  is  sufficiently  complex  to  war- 
rant a numerical  solution.  We  will  derive  now  the  general  expression 
for  the  efficiency  to  be  used  in  this  computation,  and  then  discuss 
the  results. 

Because  the  up-conversion  process  can  be  phasematched  at  certain 

, one  must  start  from 
an  expression  for  the  efficiency  thaL  Lakes  into  account  saturation. 


combinations  of  the  variables  An,  , g , and 


Such  an  expression  can  easily  be  found  by  analogy  with  the  case  of 
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crystal  up-converters,  where  a second  order  nonlinearity  is  used. 
For  metal  vapors,  the  corresponding  result  is: 


K.K  + (Ak/2)£ 
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The  expression  for  the  mismatch  Ak  is  the  same  as  for  the  active 
imaging  case. 

As  an  example,  we  will  calculate  the  response  versus  angle  for 
a cesium  up-converter  using  the  upper  level  of  the  p doublet.  In 
the  calculation,  we  will  neglect  the  contribution,  both  to  the  ef- 
ficiency and  the  k-vector  mismatch,  of  the  lower  level  of  the  p 
doublet.  Since  the  upper  level  of  the  doublet  has  two-thirds  of  the 
oscillation  strength,  this  should  result  in  only  a small  error  for 
detunings  from  the  upper  level  less  than  about  half  the  spacing  be- 


tween the  two  levels,  which  is  1 1 cm  for  the  case  of  that  transi- 
tion. So  for  bandwidths  up  to  about  cm  ^ 1 cm  ^ on  each  side 


of  the  level),  the  error  introduced  by  making  this  simplification 
should  be  small.  Because  of  absorption  at  the  sum  frequency  for  very 
small  detunings,  a small  band  of  frequencies  has  to  be  excluded  in  the 
integration  over  frequency.  For  the  case  of  the  cesium  up-converter, 
for  a density  of  1 ^ cm  ' and  a cell  length  of  1 mm.  the  width  of 
this  band  is  about  cm  ' . This  is  small  compared  to  the  total  band- 
width, and  its  exact  value  has  little  effect  on  tin  frequency  averaged 
efficiency.  We  will  also  neglect  the  effect  of  a small  population  in 
the  s level.  When  these  approximations  are  made,  the  complete  ex- 
pression for  the  mismatch,  including  the  contribution  due  to  the  angle 
of  the  IR  vector  is: 

N / 1 .it  \ 10"  \ 

2k  = ( - . l I t i — ~ tf  (3  -3^ ) 

2Cpdi  \ Aii  ' c c 

Using  the  expressions  (3-33)  and  (3-  1 > the  efficiency  can  be  calcu- 
lated for  every  set  of  the  variables  V'.  , ' . and  , and  from 

this  a frequency  integrated  response  versus  angle  1 can  be  found. 

The  expression  for  the  efficiency  involves  a very  large  number  of 
parameters,  however,  and  it  is  convenient  to  rewrite  it  in  a form 
that  only  contains  a small  number  of  key  variables.  We  will  rewrite 
the  expression  . in  terms  of  the  calculated  or  measured  effi- 
ciency in  one  coherence  length,  for  the  special  ease  where  the  sum 
frequency  comes  close  enough  to  the  upper  level  such  that  it  sets 
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Substituting  the  expression  for  in  '.33 we  finally  find: 
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where,  using  the  parameters  of  the  Cs  up-converter: 
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1 V 


( f i) 


As  is  clear  from  its  definition,  is  the  unsaturated  efficiencv  in 

one  coherence  length,  for  an  IK  frequency  such  that  tlu  ctlienno  lcr  th 
is  set  by  the  upper  level.  For  a given  . the  intensity  required  to 

obtain  that  efficiency  is  found  from  3- 

The  results  of  the  computation  for  the  cesium  up-converter  are 
shown  in  Fig.  3,  for  the  sets  of  parameters  shown  there.  Frequencies 
with  detunings  of  less  than  15  cm  * were  excluded  from  tlu  integration. 
If  the  computations  are  done  for  a different  IK  bandwidth,  one  finds 
that  the  resulting  solid  angle  is  such  that  the  product  of  bandwidth 
and  solid  angle  is  nearly  the  same.  Changing  the  number  density 
changes  this  product  proportionally.  This  proportionality  can  also  4 
be  derived  analytically,  if  some  approximations  and  assumptions  are 
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Photon  conversion  efficiency  of  cesium  up-converter  as  a function  of  IR  incidence 
angle  and  pump  phasematching  angle,  for  given  bandwidth . 


made.  If  the  angle  cp  is  chosen  such  that  only  the  term  containing 


and  the  term  containing  9 remain  (i.e.,  the  coherence  length 
is  set  by  the  upper  level),  the  efficiency  in  one  coherence  length 
is  proportional  to: 
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Clearly,  for  angles  9 and  detunings  An,  such  that  the  second  tenn 

in  the  denominator  is  less  than  one  (in  absolute  value),  the  envelope 

, 2 

of  the  efficiency  is  almost  constant  (the  sin  factor  with  average 
value  0.5  oscillates  under  this  envelope).  By  taking  for  9 and  An, 
their  maximum  values,  we  find  the  condition: 
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this  can  be  written  as 


AS2  x B ~ 2aX.  cc  N 
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The  numerical  vaLue  of  the  product  obtained  this  way  is  somewhat  low. 
The  general  behavior  predicted  by  it  is  also  observed  in  the  exact 
numerical  solution,  however. 

b . Resolution 

The  general  expression  for  the  resolution  is  derived  in  Appendix 
B.  For  the  solid  angle  in  this  expression  one  uses  the  value  found  in 
the  preceding  section.  For  the  cesium  example  this  angle  is  larger 
than  the  maximum  value  allowed  by  condition  B.  in  Appendix  B.  One 
should  remember,  however,  that,  except  for  a few  1R  frequencies  where 
the  process  is  phasematched,  the  generation  happens  in  a length  (the 
last  coherence  length)  which  is  considerably  shorter  than  the  physical 
length  of  the  medium.  For  the  same  example  as  above,  the  coherence 
length  for  large  detunings  equals  0.2  mm  [from  which  is  five 

times  shorter  than  the  length  of  the  vapor  zone.  It  is  this  length 
that  should  be  used  in  the  condition  on  the  solid  angle,  and  if  this 
is  done,  there  is  no  violation. 

c . Relation  between  efficiency,  resolution,  and  number  density 

The  efficiency  in  one  coherence  length  is  proportional  to  (since 
NLc  is  independent  of  N ) : 
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Multiplying  both  sides  by  R"~  , and  using  (3.1*2),  we  find,  for  a given 

IR  bandwidth 


(3.W) 

which  is  the  same  result  as  for  the  active  imaging  case.  It  says  that 
one  wants  to  operate  at  the  maximum  pressure  allowed  by  limiting  ef- 
fects . 

D.  PHENOMENOLOGY  OF  TWO-PHOTON-RESONANT  UP -CONVERTERS . EXAMPLE  OF 
2.9I*  p UP-CONVERSION  IN  CESIUM  VAPOR 

1 . Efficiency 

The  expressions  for  the  efficiency  of  two-photon-resonant  fre- 
quency converters  were  derived  in  Chapter  II.  They  are  repeated  here 
for  convenience: 


p O p p p 

£ X R <x  P Aft  « P N 

u 


Doppler  broadened  regime 


£ 


2 1*  2.2  2 

T]  K 'L  I 

1.8  x -s-  0 1 c x 


7 


D5 


- 102  - 


or 


or 


^0^ 

c » ) i S U 1 c. 

C - 4.4  x X sine 

V 


(?) 


(5.^5) 
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where 


K = Nm01m12^23^30 

1 bn?  (o>1Q  - aiQ)  (cUjo  - ) 

7n  = Doppler  linewidth  = F.W.H.M./2.35 


= pressure  broadened  linewidth  s F.W.H.M. 


$ = laser  bandwidth  = 1 .2  X F.W.H.M. 


(2.15) 


We  will  apply  these  expressions  to  the  case  of  .9^*  p up-conversion  in 
cesium  vapor.  Figure  9 shows  the  energy  levels  of  cesium.  The  rele- 
vant atomic  parameters  are  tabulated  in  Table  I. 
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CESIUM  ENERGY  LEVELS 


-Energy  levels  of  cesium  and  relevant 
wavelengths  for  IR  up -converter . 


TABLE  I 
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],  lengths  in  [m],  and  the  atomic  density 


Limiting  Effects 

a . Molecular  absorption  at  the  pump  frequency 

The  absorption  cross  section  of  Cs  dimers  has  been  measured  by 
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Kostin  and  Khodovoi.  Their  results  are  shown  in  Table  II.  The 

- NL 

transmission  through  a length  L equals  e . As  an  example, 


I ) 


TABLE  II 

ABSORPTION  CROSS  SECTIONS  OF  Cs  DIMERS  AT  1.079  D 




Temperature  [ C] 

N [ cm  ] 

Cs 

N [cm  ' ] 

Cs2 

■ [ cm'  ] 

320 

5.9  X 101' 

rH 

0 

T-* 

X 

i».:  X 1 _1 

380 

2.0  x 101 

— 

. x lo15 

.0  ■ 10_1 
— 

107  ' 


at  a temperature  of  JO  C,  the  transmission  for  a 2 mm  long  cell  is  9 0 

At  higher  temperatures,  the  cross  section  has  unfortunately  not  been 
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measured.  Workers  at  NRL  claim  that  the  dimer  absorption  is  reduced 
by  overheating  the  vapor  using  a side-arm  cell,  where  the  temperature 
of  the  up-converting  region  can  be  much  higher  than  the  side-arm  tem- 
perature. Unfortunately,  at  the  time  of  this  writing,  no  quantitative 

22 

data  were  available  yet  on  this  effect.  Kostin  and  Khodovoi  have 

observed  clarification  of  the  absorption  bands  of  Csn  at  1 . ■ inten- 

7 , 2 

sities  of  1 W/cm  and  up.  Unfortunately,  these  intensities  are  con- 
siderably higher  than  those  used  in  the  up-converters. 

Additional  measurements  should  be  made  on  the  absorption  cross 
section  of  Cs^  dimers  at  high  temperatures.  They  should  oe  done  in  a 
side-arm  cell  such  that  the  effect  of  overheating  the  vapor  can  be 
studied . 

Table  111  gives  the  coefficients  for  the  vapor  pressure  depen- 
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dence  on  temperature  of  Cs  and  Cs0.  ' In  Fig.  1 , the  number  density 

C. 

of  Cs  and  Cs^  are  plotted  versus  temperature  between  '•  lC  and  ■ °C . 
In  Fig.  11  the  ratio  of  the  atomic  to  the  molecular  number  density  is 
plotted  versus  temperature. 

b . Two-photon  absorption  of  the  pump  radiation 

The  expression  for  the  two-photon  absorption  rate  was  derived  in 
Chapter  II  (2.55)*  For  the  cesium  up-converter: 
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cesium  versus  temperature. 
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where  I is  in  [w/m  ] and  & is  in  [cm  ].  Using  this  expression, 
the  attenuation  of  the  pump  beam  can  be  calculated: 
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dz 
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where 


m w ■ i2 


The  solution  to  this  equation  is: 


1(0 


1 + ONLI 


(5-53 


where  I is  the  intensity  at  the  entrance  of  the  cell.  For  the  case 
of  the  cesium  up-converter,  the  transmission  of  the  cell  is: 


1(0 
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By  combining  the  expressions  for  the  efficiency  and  for  the  two- 
photon  attenuation  of  the  pump  beam,  the  maximum  efficiency  allowed  by 
this  effect  can  be  calculated.  We  will  limit  the  quantity  ONLI  to 
a value  of  0.5,  which  corresponds  to  a transmission  of  67".  For  the 
case  of  phasematched  up-conversion,  the  resulting  maximum  efficiencies 
are : 


For  the  cesium  up-converter: 
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As  an  example,  for  5 = 0.12  cm  * , a detuning  of  . cm  ' and  a pres- 
sure broadened  linewidth  of  0.1  cm  ^ , we  find  a maximum  efficiency  of 
8.3,  corresponding  to  a photon  efficiency  of  1.3.  This  means  that  a 
photon  efficiency  of  107  would  be  possible  in  this  example. 
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Saturation  of  the  two-photon  transition 


c . 


The  dipole  moment  at  the  sum  frequency  is  proportional  to  the 
population  difference  between  the  ground  level  and  the  resonant  level. 
The  intensity  thus  has  to  be  limited  such  that: 


X r 


v 0.5 
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where  r is  either  the  pulselength  or  the  lifetime  of  the  resonant 
level,  whichever  is  shorter.  For  the  cesium  up-converter,  for  5 
0.12  cm  1 , a pulselength  of  100  nsec,  and  a lifetime  of  the  's-ip 
transition  of  52  nsec,  we  find  a maximum  allowed  intensity  of  5x1 
W/cm  . 


d . Atomic  absorption  at  the  sum  frequency 

For  detunings  larger  than  a few  atomic  linewidths,  the  single 
photon  absorption  cross  section  is  given  by:^ 
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For  the  cesium  up-converter,  7 - 5.27  X 10  ' \ N [cm 
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transition  6s‘ S-7p'  P.  /n . The  transmission  of  the  cell  is 
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with  At' . in  [cm  ^ ] and  N in  [m  ]. 
Examp  le : 
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e . Atomic  absorption  at  the  IR  frequency 

Even  a small  population  in  the  resonant  level  will  result  in  a 
considerable  attenuation  of  the  IR  field,  since  this  field  is  almost 
resonant  with  the  allowed  transition  from  that  level  to  the  upper 
level.  In  Appendix  D the  pressure  broadening  of  this  transition  is 
estimated  for  the  case  of  the  cesium  up-converter.  The  result  is 
y ~ u x 10  X N [cm  1 ] , with  N in  [m  ] . Using  this,  we  find 
for  the  cesium  up-converter: 


oNL  ~ 9.0  x 10 
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where  f is  the  fraction  of  atoms  in  the  s level,  in  percent.  In 
deriving  this  result,  use  was  made  of  the  fact  that  the  broadening  is 
almost  completely  due  to  the  atoms  in  the  ground  state,  since  only  a 
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very  small  fraction  is  allowed  in  the  s level,  while  the  absorption 
is  done  by  the  small  fraction  of  atoms  that  is  in  that  level. 


Example : 
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An, 

■ 

-»  transmission 
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2 k -3 
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2 x 10_:' 

20  cm  ^ 


f . Breaking  of  phasematching  due  to  a non-zero  population  in  the 

resonant  level 

If  the  up-conversion  process  is  self-phasematched,  and  no  popula- 
tion is  present  in  the  resonant  level,  then  it  will  be  phasematched  at 
any  number  density  since  all  the  contributions  to  the  mismatch  are  pro- 
portional to  this  number  density.  This  is  not  the  case  if  the  third 
technique  for  phasematching  is  used,  since  in  this  case  the  angle 
depends  on  the  number  density  according  to  ■ I'1  Section  C,  the 

effect  of  a non-zero  population  in  the  resonant  level  on  the  mismatch 
'k  and  the  phasematching  frequency  (3*25)  and  angle  was  dis- 

cussed in  considerable  detail.  Since  the  population  of  the  resonant 
level  is  a function  of  time,  a distinction  has  to  be  made  between  the 
case  where  the  pulselength  is  shorter  than  the  decay  time  of  the  reso- 
nant level,  and  the  case  where  the  pulselength  is  much  longer  than  that 
time,  such  that  steady-state  conditions  are  present  during  most  of  the 
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pulse.  In  the  first  case,  the  population  of  the  resonant  level  has  to 
be  limited  to  a value  such  that  at  the  end  of  the  pulse  the  process  is 
still  phasematched . In  the  second  case,  a larger  population  can  be 
tolerated  because  its  effect  can  be  offset.  However,  since  this  pop- 
ulation (steady-state  value'  depends  on  the  pump  power,  the  effect  of 
pulse-to-pulse  fluctuations  of  the  laser  power  must  be  considered. 

Also,  in  the  second  case  the  beam  profile  of  the  pump  is  important, 
since  the  population  in  the  resonant  level  will  spatially  vary  ac- 
cording to  the  square  of  the  local  intensity,  such  that  the  process 
is  only  phasematched  over  part  of  the  pump  beam,  unless  beam  apodiza- 
tion  is  used.  We  will  now  study  both  these  cases,  for  the  example  of 
the  .0';  P cesium  up-converter. 

For  pulses  shorter  than  the  decay  time  of  the  resonant  level,  the 
phase  mismatch  iias  to  be  less  than  <T  at  the  end  of  the  pulse,  i.e.: 

N. . I.l£  m 

■rn  -'8-  --1  < 7T  (3 .<’3) 

2e  ctt&n 

0 3 

Substituting  the  Cs  parameters,  this  requires: 
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with  N in  [m  '} ) , L in  [m],  and  An  in  [cm  ^ ] . Since 
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With  this  fraction  of  atoms  f in  the  resonant  level,  the  process  will 
self-phasematch  at  a detuning  given  by  (3.25)  or,  for  a given  detuning, 
the  phasematching  angle  Cp  will  be  given  by  (3.28  . In  this  case,  a 


larger  fraction  can  be  allowed  than  in  the  previous  case,  resulting  in 
higher  efficiencies.  As  an  example,  when  the  fraction  f is  increased 
to  0.6  from  its  0.2  value  in  the  other  case,  the  square  of  the  intensity 
can  be  increased  by  a factor  of  three,  resulting  in  an  efficiency  three 


times  higher.  The  process  will  in  this  example  be  self-phasematched  at 
a detuning  of  10.5  cm  ^ instead  of  16.3  resulting  in  an  additional  in- 


crease in  the  efficiency  because  of  the 
case  of  angle  phasematching  with  An  -- 


l/Ai| 

20  cm 


dependence)  or, 
, the  angle  Cp 


in  the 
will 


increase 


from  1 .3°  to  2.1°. 


Because  the  steady-state  population  is 


proportional  to  the  square  of  the  intensity,  one  has  to  watch  in  this 


case  the  pulse-to-pulse  stability  of  the  pump  laser,  and  also  the  ef- 
fect of  a non-rectangular  beam  profile.  The  condition  for  the  pump 
stability  can  be  found  as  follows.  If  there  is  a change  'll  , 
caused  by  a change  AL  , the  condition  AcL  < 7T  is  the  same  as 

except  that  the  population  N,-,  is  replaced  by  the  change 
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in  that  population.  In  terms  of  -XI 


, the  condition  is: 
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Using  the  same  parameters  as  in  the  previous  example,  we  find  for  the 
sel  f-phasematched  case  (A:\.  1 .5  cm  * MI  5 . and  for  the 

angle  phasematched  case  (An.  20  cm  * .M/I  '.5  . These  sta- 

bilities are  easily  obtainable  with  present-day  solid-state  lasers. 

Since  ‘'ll  /N  M/I  , we  can  also  find  the  time  after  which  the 

s s 

process  becomes  phasematched,  using  (y.t  . For  M/l  1 1 , we 
find  that  the  mismatcii  is  less  than  n after  t 1 .t  > r,  ^ . For 
the  cesium  up-converter,  this  means  that  the  pulse  length  should  be 
considerably  longer  than  3 nsec.  From  the  value  of  \l/l  , we  can 
also  calculate  the  maximum  radius  at  which  the  process  will  still  be 
phasematched  in  the  case  of  a Gaussian  pump  beam.  As  an  example,  for 
Al/l  - 10  , the  maximum  radius  is  equal  to  . ;3  w , while  the  ef- 

O 

ficiency,  which  is  proportional  to  I'  , only  drops  to  5'  of  its 
value  in  the  center  for  r 0.'*2  w . In  such  a case  one  might 
want  to  apodize  the  beam. 


g . Breaking  of  phasematching  due  to  quadratic  Kerr  etlcct 

The  change  in  the  index  ef  refraction  at  the  pump  frequency, 
caused  by  the  quadratic  Kerr  effect,  was  calculated  in  Chapter  II. 
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It  was  found  that  the  maximum  value  of  this  change  is  given  by: 


(n  - 1 


r 1 M 


max 
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For  the  cesium  up -converter : 
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(n  - 1)  = 2.7  x 10'^  — 

max 
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where  N is  in  [m  ] . I in  [W/m‘  ],  and  y in  [cm  ].  The  phase 
mismatch  over  a distance  L is  then  given  by: 


2<n  (n  - 0L 


.'kL 


For  the  cesium  up -converter: 
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Example : 
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For  this  set  of  parameters,  we  find: 


(n  - 1)  = 5.1+  x 10‘ 

max 


AkL  ~ 0 . 1 « 7T 


Because  the  IR  field  is  very  weak  in  imaging  applications,  the 
Kerr  effect  at  the  sum  frequency  can  be  neglected.  Ths  would  not  be 
the  case,  however,  if  one  would  try  to  deplete  the  pump  beam  with  a 
high  power  IR  beam. 


h.  Thermal  defocusinc 


If  the  pump  beam  is  not  apodized,  the  population  of  the  resonant 
level,  which  has  the  same  spatial  dependence  as  the  square  of  the  in- 
tensity, will  create  a lensing  effect  for  the  incident  waves.  This 
happens  because  for  the  IR  waves  -\n  « N , while  for  the  pump  and 
the  sum  waves  .Vi  « ,\N  g - N . . When  the  pump  beam  is  split  into 

two  beams  which  then  interfere  at  a small  angle  ; , the  transversal 

variation  of  the  total  pump  field  will  be  sinusoidal,  with  period: 


\ 

X = -E- 

P CM 


The  population  of  the  resonant  level,  and  thus  the  index  of  refraction 
will  vary  sinusoidally  in  this  case,  with  a spatial  period  one-quarter 


that  of  the  pump  field: 


5. 


The  incident  waves  thus  see  a periodic  phase  grating.  The  Kerr  con- 
tribution to  the  index  of  refraction,  which  is  proportional  to  the 
intensity,  varies  with  a spatial  period: 


X 

P, 


Kerr 


X 
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creating  another  phase  grating.  We  will  now  estimate  the  importance 
of  these  effects  for  the  example  of  the  Cs  up-converter. 

The  effect  of  the  Gaussian  profile  of  the  population  in  the  reso- 
nant level  can  be  estimate d as  follows.  According  to  Kef.  1 A the  con- 
dition that  has  to  be  satisfied  to  avoid  thermal  defocusing  is: 
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w 
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where  An  is  the  change  in  the  index  of  refraction  caused  by  the  atoms 
in  the  resonant  level,  and  the  factor  comes  from  the  dependence  on  the 
square  of  the  intensity.  For  the  1R  waves,  for  the  example  of  the  cesium 
up-converter : 
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Example : 


, 2l+ 

N = 

1C  m 

f = 

0.2 

20  cm 

-»  7n  = 

7 x l 

For 

L - 0.2  cm  , K 10'  , 

the 

radius 

w is 

P 

on  the  order  of  1 cm. 

Th  i s 

yields  */  .In  X L/w 

«.  10 

_ 

« 1 

. For 

the  pump  and  sum  waves 

the  matrix  element  is  smaller  and/or  the  detuning  larger,  and  as  a re- 
sult condition  (3.75)  is  satisfied. 

The  effect  of  the  Kerr  contribution  to  the  index  of  refraction  at 
the  pump  can  also  be  estimated  using  (3-  5>  where  n is  now  due  to 
the  quadratic  Kerr  effect.  Using  (;  . . with 
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q / 2 

1 1 W 'm 

. 5 cm 


one  finds  n,  5 < 1 . Since  this  value  is  even  much  smaller 

Kerr 

• 

than  the  value  of  . \n  » the  effect  should  he  negligible. 

s 

The  problem  of  the  diffraction  of  a plane  wave  by  a sinusoidal 
phase  grating  is  considered  in  detail,  in  Ref.  . In  the  Raman -Nath 


regime,  i.e.,  for 
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L 

max 


k\ 


0.'" 


where  X is  the  period  of  the  grating  and  is  the  wavelength  of 

the  incident  light,  the  diffraction  into  higher  order  modes  is  small 

when  \kL  1 . Note  that  this  is  the  same  condition  as  tiie  one  for 

breaking  phasematching  due  to  atoms  in  the  resonant  level . Since  the 

length  L is  only  a few  microns  for  typical  values  of  , we  are 

° max 

in  the  Bragg  regime.  In  this  regime,  considerable  coupling  to  higher 
order  modes  only  exists  when  a wave  is  incident  at  the  Bragg  angle, 
i.e.: 


0 ~ 

?X 


For  the  IR  waves,  this  condition  can  be  written  as: 
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This  angle  is  much  larger  than  the  acceptance  angle  of  the  up-conver- 
ter, and  rays  incident  at  such  a large  angle  can  be  stopped  using  a 
mechanical  stop.  Also,  in  the  Bragg  regime,  the  condition  for  sig- 
nificant scattering  into  higher  order  modes  is  again  AkL  ~ 1 

3 . Optimization 

We  will  address  the  problem  of  finding  the  optimum  values  for  the 
atomic  number  density  N and  the  length  of  the  cell  L , for  a given 
efficiency  and  resolution  . Because  active  imaging  looks  much  more  prom- 
ising than  passive  imaging,  when  compared  with  competing  techniques.^ 
only  that  case  will  be  discussed.  The  procedure  will  be  explained  on 
the  hand  of  the  . p Cs  up -converter . 

In  Table  IV  the  most  important  limiting  effects  are  summarized, 
together  with  the  quantity  they  bound.  To  solve  the  problem  we  need 
in  addition  the  expressions  for  the  efficiency  [ and  .bl  and 

the  expression  for  the  number  of  resolution  elements  ‘ . • . Because 

at  high  number  densities  the  atomic  Linewidth  becomes  pressure  broad- 
ened, such  that  the  dependence  of  the  efficiency  on  number  density  is 
linear  instead  of  quadratic,  the  Doppler  and  the  pressure  broadened 
case  have  to  be  considered  separately.  The  optimization  procedure  is 
very  similar,  however,  for  both  cases.  By  limiting  the  molecular  ab- 
sorption and  the  population  in  the  resonant  level  *,to  avoid  breaking 
phasematching'  a maximum  value  can  be  found  for  the  quantity  NP  K 
By  substituting  this  value  in  the  expression  for  the  efficiency,  a 
solution  to  the  problem  is  obtained  for  both  cases.  This  solution  is 
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then  substituted  in  the  remaining  conditions  of  Table  IV  to  check 
whether  they  are  not  violated.  In  view  of  the  results  on  thermal  de- 
focusing,  in  applying  the  condition  on  breaking  of  phasematching,  we 
will  limit  the  population  in  the  resonant  level  to  a value  such  that 
phasematching  is  not  broken  when  the  IR  frequency  is  tuned  to  the 
phasematching  value  for  f , or,  for  given  IR  frequency,  when 

the  angle  is  adjusted  to  its  value  for  f . In  the  self- 

phasematched  case  , a higher  population  can  be  tolerated 

in  the  resonant  level  ii  long  pulses  are  used,  but  the  power  required 
to  obtain  a given  efficiency  will  be  found  to  vary  only  as  the  square 
root  of  the  increase  in  population,  making  it  a small  effect. 

To  limit  the  molecular  absorption,  we  require: 


which  corresponds  to  transmission  of  the  pump.  Since  we  are  not 
interested  in  the  molecular  number  density,  we  replace  it  by  the  prod- 
uct of  tin.  atomic  number  density  and  some  proportion.!  1 itv  constant, 
which  can  be  done  over  a small  temperature  range  Fig.  11  . Around 
C.  this  ratio  is  about  . The  condition  . can  then  be 
written  as: 


NL 
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To  avoid  breaking  of  phasematching,  we  require  (5.65): 


5 X 10 


30 


(3.32) 


By  combining  this  with  the  expression  for  the  resolution  (3.29: j this 
condition  can  also  be  written  as: 


NP2  . 

— < 1.75x10-  2 (3.:3) 

LR"  r 


Finally,  the  combination  of  this  result  with  the  condition  on  molecu- 
lar absorption  yields  the  important  relation: 

NP  / B&D,  \ 

— ^ 1.5  X 1 L I 2 (3-3*0 

R \ CTT  / 

We  now  have  to  consider  separately  the  cases  of  Doppler  broadening 
and  pressure  broadening.  The  Doppler  width  (F.W.H.M.  for  Cs  is 
cm  ^ at  a temperature  of  UUo°C.  In  Appendix  D the  pressure  broadening 
of  the  linewidth  of  the  os-'s  transition  is  estimated  to  be  cm  ^ 

at  a number  density  of  6 X 10^  cm  . Thus  for  densities  below  1 ^ 

cm  J we  are  in  the  Doppler  broadened  regime,  while  for  densities  higher 
than  this  value  we  are  in  the  pressure  broadened  regime. 

a.  Doppler  broadened  regime:  N 1 • 1 * cm* 

Combination  of  the  expression  for  the  efficiency 
expression  for  the  resolution  (3.291  yields: 
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and  the 


The  second  factor  in  this  expression  is  exactly  the  square  of  the  quan- 
tity for  which  we  already  found  an  upper  bound  •'  . . Substituting, 

we  find  the  very  important  result: 


1 . 1 


■ An, -t 
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It  says  that  in  the  Doppler  broadened  regime,  there  is  a maximum  ob- 
tainable efficiency,  which  is  independent  of  the  resolution.  Hie  com- 
bination of  this  result  with  (p  . . . 1 , and  . solves  the 

optimization  problem.  From  (p  . M it  is  clear  that  for  a given  reso- 
lution, we  want  to  maximize  N . Thus  we  want  to  use  the  maximum 
density  allowed  by  the  requirement  that  we  are  in  the  Doppler  broad- 
ened regime.  Using  N N ^ .5  1 cm  see  Fig.  . the 

relation  (p . between  power  and  resolution  can  be  written  as: 


Examp  1 e : 


?D 

20  cm 

0.013  cm  ^ 

rr 

, -21  2 

( X 10  m 

T 

50  nsec 

p 

, -1 
. 1 cm 

With 

these  parameters,  we 

find: 

'"d 

1 . •< 

p 

5 •■  ■ ■ R W 

L 

.1  cm 

A 

1.  ■ 1 " 

R cm 

The 

remaining  limiting  effects  can  now 

be  checked. 

Since  1 

are 

fixed,  we  find,  independent  of  the 

resolution 

Table  IV 

Absorption  at 

s • 1 : 

NL 

. 1 

Absorption  at 

ir  • : 

NL 

Absorption  at 

P • : 

iINL 

. 

b. 

Pressure  broadened  regime:  N ^ 

1 cm 

Combination  of  the  expressions  for 

tlie  efficiency  A 

resolution  . 1 yields: 


- 


Since  7 > 7^  , and  NP/R  is  a constant,  we  can  say  immediately  that 

the  efficiency  in  the  pressure  broadened  rerime  will  be  less  than  in 


the  Doppler  broadened  regime.  This  does  not  imply,  however,  that  the 
maximum  pressure  allowed  by  the  condition  7 < 7 is  always  the  opti- 
mum operating  condition.  If  the  desired  efficiency  is  less  than  the 
maximum  value  allowed  in  the  Doppler  regime,  then  the  required  power 
in  the  pressure  broadened  regime  will  go  down  proportional ly  to  the  ef- 
ficiency decrease,  since  £ « l/n  and  P - r/n  , while  in  the  Doppler 
regime  it  will  only  go  down  as  the  square  root  of  this  decrease  . 5 • 
We  will  now  derive  a relation  between  number  density  and  efficiency, 
which  is  valid  in  the  pressure  broadened  regime  but  can  only  be  used 
for  efficiencies  less  than  the  maximum  efficiency  obtainable  in  the 
Doppler  regime.  In  Appendix  D,  the  linewidth  is  estimated  to  be: 


~ 5 X IQ-"''  N [cm'1  ; 


where  N is  in  [m] . Substitution  in  . yields: 


l.t  X 1 


XX.-'N 


(t) 


After  combining  this  with  . , we  find  the  desired  result: 


5.6  x 10  l 


Example : 


result 

with  5 . 

problem 

• 

II 

<1 

20  cm 

a = 

7 x 10" 

T = 

r'  nsec 

£ = 

. If  cm 

e = 

1 . 

R = 1 


With  these  parameters,  we  find: 

N = 

P = 

L 

A - 6. 


5 X 101  cm" 


2.5  '■  10"  w 

1 .5(  cm 

A 

cm' 


3 • •! 

. 

. 1 


The  remaining  limiting  effects  can  now  be  checked: 


Absorption 
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.<1  : 

NL 

Absorption 
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ir 

. : 

•NL 

. 5 

Absorption 
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'0 

P 

C(N  LI 

The  number  density  of 

5 X 

i 17 

cm  cor res 

ponds  to  c 

i temperature  ol 

about  y.  C,  where  the  molecular  cross  section  is  not  exactly  known 
This  density  is  only  a factor  of  . higher  though  than  the  maximum 


1 


i 


Combining  ,3 •T’}  and  (3.78): 


7 - 9.8  > l 5 x 


A x P 


M ■ 


L x P 


■ 


Because  of  the  L'  /A  dependence  of  the  efficiency,  it  is  optimum  to 
use  as  long  a length  as  possible  for  low  resolutions  . As  the  resolu- 
tion is  increased,  however,  the  area  A will  also  increase,  and  at 
some  "critical"  resolution  the  area  becomes  equal  to  some  maximum 

value  A . If  the  resolution  is  increased  beyond  that,  the  length 

max 

has  to  be  decreased,  and  formula  . 7 must  be  used  with  A A 

max 

Clearly  the  critical  resolution  is  given  by: 


A 

„ , , b max 

2.3  x 10  x 


We  will  now  make  a general  comparison  between  vapor  and  crystal 
up-converters.  Two  separate  cases  must  be  considered,  depending  on 


whether  R R • We  will  start  by  comparing  the  up-converters  at 

the  maximum  efficiency  • .86  of  the  vapor  up-converter.  After- 
wards we  will  make  some  comments  for  the  case  where  L'  ''  . For 

the  cesium  up-converter,  the  same  parameters  are  used  as  in  the  section 
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We  will  not  discuss  this  problem  here  because  it  is  of  little  interest 
for  imaging  applications  where  the  highest  possible  efficiency  is  needed 


because  of  the  small  signals  received.  The  design  parameters  for  a 
crystal  and  a vapor  up-converter,  with  £ - £ and  R = 10^  , are 

summarized  in  Table  V. 

E.  EXPERIMENTAL  RESULTS 
1 . Up-Conversion  of  : ,9--  n Images 

a . Discussion  of  the  experimental  set-up 

The  set-up  for  the  imaging  experiment  is  shown  in  Fig.  12.  The 
Nd:BeL  rod  used  in  the  Q-switched  laser  has  a line  which  comes  within 
about  1 X of  the  two-photon-resonance,  and  it  could  easily  be  tuned  on 
resonance  with  a Fabry-Perot  etalon  inside  the  cavity.  The  beam  coming 
out  of  the  laser  was  split  into  two  parts.  About  * of  the  power  was 
focused  into  a CDA  doubler,  the  output  of  which  was  used  to  drive  a 
LiNbCL.  parametric  oscillator.  The  pump  and  the  signal  of  the  0P0  were 
filtered  out  with  a germanium  filter,  while  the  idler  wave  was  sent 
through  a spatial  filter  and  then  collimated  to  a diameter  of  about 
3 cm.  The  plane  wave  idler  radiation  was  then  sent  through  a quartz 
slide  on  which  the  object  was  either  taped  (Figs.  13  and  l^t  or  painted 
(Fig.  15)-  The  up-conversion  was  done  using  a Type  III  optical  system. 
The  object  and  the  center  of  the  Cell  were  in  the  two  focal  planes  of  a 
20  cm  quartz  lens,  while  the  cell  and  the  film  were  in  the  two  focal 
planes  of  a V:  cm  quartz  lens.  The  magnification  was  thus  M \ 

X 50/2 ' ~ . The  remaining  ; t"  of  the  pump  was  blown  up  to  an  area 
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DESIGN  EXAMPLES  OF  A CESIUM  AND  A LiNbO., 
UP-CONVERTER  FOR  2. M RADIATION 
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Experimental  set-up  for  IR  up-conversion  in  cesium 


The  pump 


of  about  0.1  cm'  ~,  and  used  as  the  pump  of  the  up-converter, 
and  IR  beams  were  combined  using  a high  reflector  for  1 .06  p radiation. 
High  speed  polaroid  film  (3000  ASA)  was  used  to  record  the  images.  The 
shutter  of  the  camera  was  opened  for  1 / 1 2 5 sec,  and  provided  a short 
circuit  that  was  used  to  trigger  the  laser.  A blue  filter  was  put  in 
front  of  the  camera  to  black  out  most  of  the  light  of  the  laser  flash- 
lamp  . 

b . The  images 

The  cell  consisted  of  two  sapphire  windows  (manufactured  by  EIMAC), 
heli-arced  inside  two  Varian  flanges  with  a distance  between  them  of 
about  2 mm  i,Fig.  1>  '.  The  corresponding  solid  angle  given  bv  ttX. 

is  equal  to  k.C  x 10  ' steradians,  and  the  area  of  the  object  "seen" 

2 

by  the  up-converter,  given  by  AI  x f , is  equal  to  1.0  cm'  . The 

u 

2 

area  of  the  pump  beam  was  0.1  cm  , and  thus  the  theoretical  resolution 
was : 


A /2  X Aft 

R -2 ~ 2670 

X1 

1 

Figures  13  to  15  show  photographs  of  three  objects  and  the  correspond- 
ing images.  Note  that  all  the  images  are  clearly  resolved,  also  the 
second  one,  which  shows  about  1000  rectangles.  The  uniformity  of  the 
images  is  not  very  good,  but  this  is  most  likely  related  to  imperfect 
filtering  of  the  illuminating  IR  beam,  and  to  imperfections  in  the 
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g.  It  --Cell  used  for  image  up-conversion. 


optics  (filters,  etc.).  We  were  not  able  to  measure  the  efficiency 
accurately  during  this  experiment  (see  the  next  section!,  but  it  was 
roughly  20/0  (power),  for  a pump  power  of  about  10  kW,  and  a number 
density  of  l.U  x 10  1 cm  J . Figure  17  shows  a photograph  of  the  cell. 


2.  Efficiency  Measurements 


Because  of  problems  with  the  imaging  cell,  the  efficiency  mea- 
surements were  done  in  a small  heat  pipe.  Since  the  cell  was  not  used 
in  the  heat  pipe  mode,  however,  there  is  some  uncertainty  as  to  what 
N and  L to  use  in  the  theoretical  expression  for  the  efficiency.  , 
Ideally,  one  would  like  to  measure  the  quantity  ,rNdz  , since  the 
efficiency  is  proportional  to  the  square  of  this  quantity.  What  one 

p 

can  measure  instead  is  /N^'dz  . This  can  be  done  by  shining  a white 
light  source  through  the  cell  and  measuring  the  width  of  one  of  the 
lines  of  the  6p  doublet.  For  a Lorentzian  cross  section,  one  can 
easily  show: 


(F.W.H.M.)' 


1 ■ I " 1 f 

x x I yNdz 


1 i 

For  the  upper  line  of  the  6p  doublet,  7 - 2.76  X 10  ^ X N , and: 


/rdz  = 5.6  x lO59  X (F.W.H.M.r 


where  the  F.W.H.M.  is  in  (cm  ].  If  a rectangular  profile  is  assumed 
for  the  number  density,  with  a length  equal  to  the  length  of  the  heater. 
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the  number  density  found  this  way  is  consistent  with  the  one  found  from 
the  temperature  of  the  outside  wall,  if  a temperature  drop  is  assumed 
across  the  wall  of  a little  less  than  ten  degrees.  The  zone  is,  of 
course,  not  exactly  rectangular,  and  probably  somewhat  shorter  than 
the  length  of  the  heater.  The  value  of  (; Ndz  ‘ was  then  approximated 
by  Lx  /N  dz  . The  wavelength  of  the  up-converted  radiation  was  mea- 
sured with  aim  Spex  spectrometer.  A detuning  of  about  9 cm  ^ from 
2 

the  p P level  was  found.  Theoretically,  with  no  atoms  in  the  s 

level,  it  should  have  been  16.3  cm  ^ . If  one  calculates  the  fraction 

f of  atoms  in  the  7s  level,  one  finds  an  average  value  f ~ .5--  , 

av 

while  at  the  end  of  the  pulse  f ~ 1.0  . Thus  the  measured  detuning 

is  in  good  agreement  with  the  predicted  value,  found  from 


An  - I6.3  - 9.7  X f [cm-1 ] (3.  5 

' Ph .match 

An  additional  fact  to  remember  in  relation  to  this  detuning  is  that  the 
IK  wavelength  was  selected  for  maximum  blue.  As  the  IR  moves  away  from 
the  phasematching  value,  closer  to  resonance,  the  efficiency  will  drop 

because  of  the  sine  dependence,  but  it  will  increase  because  of  the 

2 

Ai‘.  factor.  One  can  check  numerically,  however,  that  this  shift  is 
much  smaller  than  the  observed  one. 

The  efficiency  was  measured  using  a pyroelectric  detector.  To  mea- 
sure the  IR,  the  pump  beam  was  blacked  before  entering  the  cell,  and  no 
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TABLE  VI 


EXPERIMENTAL  PARAMETERS  FOR  2.9k  ,1  UP-CONVERSION  IN  CESIUM 


Nd:BeL  pump 

Wavelength  vacuum 

i.  790  p 

Power 

8 kW 

Pulselength 

100  nsec 

Area 

. 2 
3.1  cm 

Intensi ty 

1 W cm 

Infrared  radiation 

Wavelength 

Area 

0.1  cm‘ 

Up-Converted  radiation 

Wavelength  (vacuum 

-55  O & 

Detuning  trom  p P. 

-1 

9 cm 

Cell  parameters 

Temperature 

290°C 

Number  density 

2.0  X 1 11  cm" 

Length  of  zone 
j N‘  dz  (absorption 

2 cm 

-r. 

9..'  y 1 ^ m 

Measured 

power 

3*  > \> 

Calculated 

power 

25 

- l’i  - 

H J 


filters  had  to  be  used.  To  measure  the  generated  blue,  filters  were 
inserted  to  remove  the  transmitted  pump  and  the  IR  H.R.  at  I.06  p, 
I-58  and  5"6l  Corning  filters  . 

Table  VI  shows  a set  of  experimental  parameters,  together  with 
the  corresponding  measured  and  calculated  efficiency.  The  agreement 
between  theory  and  experiment  can  be  considered  to  be  very  good.  Part 
of  the  discrepancy  could  possibly  be  explained  by  an  imperfect  overlap 
of  the  pump  and  the  IR  beams,  but  even  without  that,  the  number  of 
theoretical  and  experimental  ingredients  of  this  experiment  was  large 
enough  to  make  such  a factor  acceptable. 
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CHAPTER  IV 


HARMONIC  GENERATION  AT  HIGH  FIELD  STRENGTHS. 

FREQUENCY  SHIFTS  AND  SATURATION  PHENOMENA* 

A.  INTRODUCTION 

26-20 

It  is  known  that  at  high  field  strengths  the  resonances  of 

interactions  between  atoms  and  electromagnetic  fields  show  intensity 
dependent  shifts.  In  two  recent  papers'  ’ ^ Harris  discussed  harmonic 
generation  and  mixing  processes  which  may  be  used  for  the  generation 
of  coherent  radiation  in  the  100  $ region.  These  processes  require 
quite  high  intensities  and  in  this  chapter  we  discuss  how  his  results 
are  modified  when  frequency  shifts  are  taken  into  account.  We  proceed 
by  solving  the  density  matrix  equations  for  a collection  of  atoms  in- 
teracting with  two  electromagnetic  fields:  a strong  incident  pump 

field  at  frequency  to  and  a much  weaker  generated  field  at  frequency 
no>  (in  what  follows  we  limit  ourselves  to  harmonic  generation,  exten- 
sion of  the  results  to  mixing  processes  is  straightforward).  The  rea- 
son why  frequency  shifts  will  appear  in  our  results  is  that  we  do  not 


Reproduced  from  Ref.  36. 
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solve  the  equations  independently,  as  is  usually  done  in  the  calcula- 
tion of  nonlinear  susceptibilities  (e.g..  Ref.  32),  but  consider  feed- 
back of  the  higher  orders  in  the  perturbation  sequence  to  the  lower 
orders.  We  demonstrate  the  technique  with  the  example  of  third  har- 
monic generation  in  a four-level  system.  Generalized  closed  form  ex- 
pressions for  generated  dipole  moment,  absorption  probabilities  and 
coherence  length  are  also  presented. 


B.  DERIVATION  OF  FREQUENCY  SHIFTS 


The  off-diagonal  elements  of  the  density  matrix  between  the  ground 
level  and  some  other  level,  for  a four-level  system,  satisfy  the  fol- 
lowing equations  (the  levels  are  labeled  from  0 to  3): 


01 


lP10P0l  + n (P00H01  + P02H21^  t 


01 


01 


02 


ia,onPno  + - (PaXo  + PaXJ  * — 


20  02  ' ^01“12  ' K03“32/ 


02 


03 


1 pox 

i®inPn-.  + ~ (0n,H'  + P_-H'  ) - — 2 


30  03  ■ R '"02“23  ■ ^00“03  T 


03 


(U.i) 


where 


"10  ■ (Ei  - Eo>/« 


H'j  = - p X (electric  field) 


= atomic  dephasing  time  for  the  0 i transition 
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We  assume  that  the  fields  are  limited  to  a value  such  that  the  change 
of  the  population  of  the  ground  level  due  to  absorption  to  other  levels 
during  the  time  the  fields  are  on  is  small,  i.e.,  p is  essentially 
constant  (this  assumption  is  discussed  in  Section  C for  various  kinds 
of  absorption  processes).  Other  on-diagonal  elements  and  off-diagonal 
elements  between  levels  neither  of  which  is  the  ground  level  are  left 
out  because  for  the  processes  considered  they  are  usually  far  from 
resonance  and  hence  have  little  effect.  The  fields  are  taken  to  be 
sinusoidal : 

i®t 

e 

E(t;u>)  = E f c.c. 

2 

ijoit 

e 

E'(t;3m)  = E'  + c.c.  (fc  .2) 

2 

The  solution  to  the  set  of  equations  (^.1)  can  then  be  written  as: 

, , \ iko>t  , . „ , , , . 

ook(t;to)  " pok  e k = 1,2,5  ^.3) 

In  the  transient  regime  the  are  functions  of  time,  slowly  varying 

compared  to  the  period  of  the  applied  field.  We  are  interested  however 
in  interaction  times  which  are  long  compared  to  the  decay  times  of  these 
transients.  If,  as  is  explained  in  Section  C,  the  power  density  is  lim- 
ited such  that  the  change  of  the  population  of  the  ground  state  due  to 
absorption  processes  is  small  during  the  interaction  time,  then  the  set 
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of  equations  (U.l)  can  be  solved  treating  p^.N  as  a constant  (quasi- 
static approach) . The  quantities  will  however  be  modulated  by 
a very  slowly  varying  envelope  due  to  the  changing  occupation  of  the 
ground  state.  Substitution  of  (U.2)  and  (1+.3)  into  (U.l)  yields,  by 
harmonic  balance,  the  following  system  of  linear  algebraic  equations 
for  the  case  where  the  generated  field  is  much  smaller  than  the  pump 
field: 


(Adx  + 

(*“2  + 
(^5  + 


‘^Ol 

l/T02 


) P, 


01 


) p, 


02 


) P, 


03 


a21P02  + a0lP00 


a!2p0l  + a32p03 


a23p02  + a03p00 


where 


An.  = cn  - ico 

j j0  J 

aij  = PijE/2ft  M = i>2,3 

aii  - ‘'uE'/2» 


The  solution  of  these  equations  is  discussed  in  Section  C.  The  reso- 
nances of  this  third-harmonic  process  are  defined  by  setting  the  real 
part  of  the  determinant  of  the  set  of  equations  (1.1)  equal  to  zero, 
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The  zeros  of  this  polynomial  in  CD  are  intensity  dependent.  For  small 

intensities  D reduces  to  D = An. Ad,  , which  is  the  well-known 
r r 1 2 3 

form  of  the  denominator  of  the  third  order  nonlinear  susceptibility. 

The  intensity  dependence  of  the  widths  of  the  resonances  is  found  from 
the  imaginary  part  of  the  determinant,  i.e.,  , where 


2avAn  to  An  - |a  | 

D.  = — i z £ ; — 

1 T T 

02  03 


Although  the  intensity  dependence  of  the  widths  of  the  resonances  of 
the  physical  quantities  studied  in  Section  C could  be  calculated  from 
the  set  of  equations  we  will  not  study  it  in  detail  in  this 

chapter.  This  dependence  was  taken  into  account  however  in  the  ex- 
ample discussed  in  Section  D. 


C.  CONVERSION  EFFICIENCY  OF  HARMONIC  GENERATION  PROCESSES 
1 

IN  THE  PRESENCE  OF  FREQUENCY  SHIFTS 


In  earlier  work'  Harris  has  derived  an  expression  for  the  con- 
version efficiency  of  a n*"''  harmonic  generation  process  where  the 
generated  frequency  comes  close  to  a single  photon  allowed  transition, 


He  shows  that  if  the  coherence  length  is  dominated  by  the  nearby 
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transition  and  if  the  pump  power  is  limited  by  n-photon  absorption  to 
that  level,  the  conversion  efficiency  is  given  by: 


where 


2find>  1 


£ = 


o.  J . /A 
On  1' 


a 


On 


J.  /a 
1 


the  single-photon  cross  section  for  absorption  to  the 
n^1  atomic  level  (this  level  is  by  definition  the 
level  closest  to  the  generated  frequency) 
the  energy  density  in  the  incident  pulse 


To  obtain  this  result  the  pump  power  density  is  chosen  such  that  the 
product  of  the  absorption  probability  per  second  and  the  pulselength 
is  0.5,  in  which  case  about  70$  of  the  atoms  are  still  in  the  ground 
state  at  the  end  of  the  interaction  (this  limit  is  somewhat  arbitrary, 
it  is  used  here  to  be  consistent  with  Refs.  JO  and  31). 

Using  the  solution  of  the  system  of  equations  (4.U),  we  find  the 
following  expressions  for  generated  dipole  moment  at  frequency  3(n  , 

single-photon  absorption  probability  per  second  to  level  3,  three- 
photon  absorption  probability  per  second  to  level  3,  and  coherence 
length: 

^aolul^a?Vl-,o 

(p(t;>f))  = — (D_  cos  3n't  + D sin  3«'t)  (k.7) 

+ D*r 

r l 
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where 


L 

c 


h 

^oI^03,2n 


D 

r 


(**.8) 


(**-9) 


k.l 


N = the  density  of  atoms 
Tl0  = 377  ohm 


For  small  field  strengths  these  expressions  reduce  to  those  cal- 
culated without  considering  feedback  effects.  Substitution  of  (U.7 
to  (!-.10)  in  the  expression  for  the  power  density  generated  in  one 
coherence  length, 

p(3“0  I M3®))  I2  (NLJ2 

= 7 — (**-11) 

A 27T 

together  with  the  requirement  W ' x (pulselength)  = 0.5  , yields 
an  expression  for  the  efficiency  of  the  process  which  is  the  same  as 
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in  the  small  field  case,  i.e.. 


2ft  • 300  1 

e = • 

°03  Ji/A 

where  it  is  now  defined,  however,  in  terms  of  a modified  detuning 
from  level  3: 


where 


A'ux. 


Ad 


1Q23' 

Ad.  Ad 

1 2 


Ad. 


a 


12 1 


D 

r 


Ad.  A) 

1 2 


a 


12 1 


2 


So  that  the  earlier  result  is  correct  if  one  replaces  the  detuning 
from  the  third  level  in  the  expression  for  the  single  photon  cross 
section  by  a modified  detuning  which  depends  on  the  intensity  of  the 
applied  field. 

31  t 

Harris  has  also  shown  that  the  conversion  efficiency  of  a n 

harmonic  generation  process  (or  n-photon  mixing ),  where  n - 1 pho- 
tons are  resonant  with  a single-photon  non-allowed  transition  to  the 
ground  state,  such  that  the  power  density  of  the  pump  is  limited  by 
(n  - 1)  photon  absorption  to  that  level,  is  given  by: 


F 


n-1 . n 1 


(l-13) 
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where  is  either  the  pulselength  or  the  lifetime  of  the  absorbing 

level,  whichever  is  shorter. 

Using  the  solution  of  the  set  of  equations  and  taking  into 

account  the  fact  that  the  sum  of  two  pump  photons  can  come  close  to 
level  2 such  that  the  linewidth  2/l\0  can  no  longer  be  neglected, 
the  two-photon  absorption  probability  per  second  is  given  by 


r(2)  = — 


02 


a 


01 1 


a 


12' 


p 

A 

2 2 

D + D7 
r i 


Substitution  of  Eqs  . (1+.7),  (^.lM,  and  (h.lo)  into  Eq , 

(2) 

gether  with  the  requirement  W'  X = 0.5  , yields 

expression  for  the  efficiency: 


(U.lM 

(it. 11),  to- 
the  following 


e = ^ L 

T1  '“OJ1  \ 

For  small  field  strengths  the  factor  in  brackets  reduces  to  1 . If  one 
limits  the  shifts  to  values  which  are  small  compared  to  the  atomic 
level  spacings,  the  allowed  power  densities  are  usually  smaller  than 
those  found  from  the  limitation  on  multiphoton  absorption.  At  such 
lower  power  densities  the  factor  in  brackets  reduces  to  1,  but  the 
efficienty  of  the  process  is  reduced  by  another  factor  which  is  given 
by: 


l°2,r 


2^3  - la12l  ^3 


.15) 


( power  density 


power  density  applied 


as  limited  by  (n-l)  photon 


n-1 


absorption 
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D . EXAMPLE 


As  an  example  we  consider  a five-photon  mixing  process  in  Li  : 


x 816  X + 67^5  X 198  X 


2 2 

The  path  used  is:  Is  — » ls2p  — » Is  ->  ls2p  — » ls2s  — > ls2p  . Level 
positions  and  oscillator  strengths  were  taken  respectively  from  Refs. 

33  and  3k . The  pump  frequencies  were  chosen  such  that  the  detuning 
from  the  unshifted  four-photon  resonance  was  1250  cm  ^ (below  the 
resonance)  while  the  detuning  from  the  unshifted  five-photon  reso- 
nance was  3750  CIn  ^ (above  the  resonance).  The  laser  pulselength  was 
assumed  to  be  30  picoseconds  and  the  linewidth  1 cm  * . Figure  18  shows 

lk  / 2 

how  at  a power  density  of  2 X 10  W/cm  the  system  becomes  four-photon 
resonant  and  also  how,  at  higher  power  densities,  the  efficiency  satu- 
rates; here  the  efficiency  was  defined  consistent  with  Ref.  3l  as  the 
ratio  of  the  generated  power  and  the  power  at  67^5  8.  Figure  18  was 
computer  generated  and  takes  into  account  the  intensity  dependent 
widths.  At  the  high  intensities  involved,  five-photon  ionization  by 
the  pump  field  at  816  X has  to  be  considered.  From, Refs.  36  and  37 
an  approximate  value  for  the  ionization  probability  can  be  calculated. 
If  the  expressions  in  Refs.  J>6  and  37  are  adjusted  for  shifts,  a lim- 

1 k / 

iting  intensity  of  about  2 X 10  W/cm  is  found,  i.e.,  if  a pulse 
with  this  intensity  is  applied,  then  about  50'  of  the  Li  ions  will 


be  ionized  during  the  pulse,  and  also  the  broadening  of  the  four-photon 
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radiation  using  a five-photon  mixing  process. 


resonance  due  to  ionization  will  be  comparable  to  the  laser  linewidth. 


It  should  be  noted  that  if  shifts  were  not  taken  into  account,  the 
limiting  power  density  (as  calculated  by  the  formulae  of  Bebb  and 

Gold  in  Ref.  36)  would  be  about  a factor  of  17  higher.  At  the  infen- 

l1*  / 2 

sity  of  2 X 10  W/cm  the  four-photon  absorption  is  well  below  the 
maximum  allowed  value.  To  keep  the  coherence  length  short  about  1 cm 
at  a pressure  of  a few  torrs),  the  power  density  of  the  c7^5  X field 
has  to  be  limited  to  about  one  part  in  ten  thousand  of  the  density  of 
the  8l6  X field.  It  is  also  indicated  in  Fig.  18  that  in  certain  re- 
gions of  the  intensity  the  dispersion  of  the  medium  can  change  sign. 

This  is  important  because  harmonic  generation  with  a focused  Gaussian 

35 

beam  in  an  infinite  positively  dispersive  medium  is  impossible.  It 
should  also  be  noted  that  the  width  of  the  resonance  in  Fig.  1“  is 
quite  small,  and  hence  a stable  laser  is  required,  typically  around  1'. 

We  conclude  that,  taking  into  account  frequency  shifts,  harmonic 
generation  is  a promising  means  of  generating  wavelengths  in  the  VUV 
and  soft  x-ray  spectral  region.  The  most  important  consequences  of 
frequency  shifts  are  that  the  pump  laser  has  to  be  tuned  to  a wave- 
length different  from  that  required  t > bring  the  system  into  resonance 
at  small  field  strengths,  that  the  conversion  efficiency  may  saturate, 
and  that  the  dispersion  of  the  medium  may  change  sign. 

E.  GENERALIZATION  TO  nth  ORDER  PROCESSES 

Using  the  technique  explained  in  Section  B it  can  be  shown  that 
til 

for  a nL  harmonic  process  (or  n-photon  mixing  , the  generated  dipole 


1 1 


moment  at  frequency  ntn  , the  single-photon  absorption  probability 

til 

per  second  to  the  n level,  the  multiphoton  absorption  probabil- 
ities per  second  to  the  (n-l)1"^  and  n^  level,  and  the  coherence 
length,  are  given  by: 
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The  efficiencies  for  the  processes  discussed  in  Section  C are: 


Case  A 
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where  o\  is  defined  in  terms  of  the  single-photon  absorption  prob- 
On 

ability  per  second  to  level  n 


Case  B: 


31 


T0n-1  K-lJ£  (,  ' °“>)l  I 

= ^ V 1 ” / 


“1  IK0,n' 
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In  these  expressions  is  defined  by: 


Dk  - “kO  - k,“ 


D1  - “10  ■ “ 


Lk-1 


k / 1 


Note:  in  all  these  expressions  the  linewidth  was  neglected. 


APPENDIX  A 


PROGRAMS  SEATPR  AND  SEATPR: 


These  programs  calculate  respectively  the  efficiency  of  tv '-photon- 
resonant  frequency  converters  versus  detuning  from  the  resonance,  and 
the  dependence  of  the  on-resonance  efficiency  on  the  pressure  broadened 
linewidth.  They  perform  a numerical  evaluation  of  . r , within  an 
arbitrary  constant.  As  they  are,  they  assume  a Gaussian  power  spec- 
trum for  the  laser,  but  this  could  easily  be  generalized  to  any  spec- 
trum, at  a point  in  the  programs  that  is  clearly  indicated.  The  input 
data  are: 


For  SEATPR 

— laser  linewidth  F.W.H.M. 

— mode  spacing 

— Doppler  width  F.W.H.M. 

— atomic  linewidth  F.W.H.M. 

— maximum  detuning 

cU 

— increment  of  (fu,. 


It  h - 


For  SEATPR. 


— laser  linewidth  (F.W.H.M. 

— mode  spacing 

— Doppler  width  F.W.H.M. 

— maximum  atomic  linewidth  F.W.H.M. 

— detuning  - 2'n.'. 

— increment  of  the  atomic  linewidth 

The  programs  are  listed  on  the  following  pages. 
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/ / S i: A T I* H -TOR  ' N 1 50 , 240 , 1 . , 1 ' , ' STAFI’A  r.RTS 

***si:rvicc  c. l a s s =q  rxr.os 

//  KXKC  WATFIV 
//(  ().  SYS  I\’  !>D  * 


SI'!'.  ItOl'T  INKS  l!  KOI’  I I'  F.D 


END  OF  LOO  I 


FAIL  W X E R R ( X 


W I'  I T I ( I UNIT  , 501) 


OOL'RLE  PRECISION  FUNCTION  HREAE(Z) 


AM.  THE  I)  A T ^ SHOULD  RE  ENTER  E 1)  IN  NAVE  MUM  R F,  K S . 

[ : U - 1 ,MO  -SPACING, 

hopple  pu  i dt  i ( fuiim)  , atom  ic-l  i nkwi  i>tm  ( fwiim  ) , 

THE  MAX . DF.TUN I NC  FROM  THE  RESON.  OM02-2*OMO, 
THE  INCREMENT  OF  THE  FREOUE'ICY  OMO. 
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CM.cn. ATim:  or  the  cuKVoLUTior;  or  the  pouek  spectrum 

FOR  C I VI'.  N v. 
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APPENDIX  B 

RESOLUTION  OF  A TYPE  III  UP-CONVERTER 


I 

[ 


We  calculate  the  magnification  of  an  up-converter  using  a Type  III 
optical  system,  and  the  effect  of  its  finite  length  on  its  resolution. 
Figure  1 shows  a diagram  of  the  complete  optical  system.  The  IK  field 
in  the  up-converter  is  the  Fourier  transform  of  the  field  in  the  object 
plane.  The  nonlinear  polarization  in  the  medium  varies  spatially  as 
the  product  of  this  IR  field  and  the  square  of  the  pump  field.  Assum- 
ing a Gaussian  pump  profile: 


dP 


CC 


u. 

lr 


2/2 


-2r‘/w 


X e 


do 


(B.r 


where  w is  the  radius  of  the  pump  beam,  and  dP  is  the  polarization 
set  up  in  a thin  slice  of  medium  with  thickness  da  . At  first  sight 
we  have  a problem  in  applying  the  results  of  Fourier  optics'  to  this 
case,  since  we  are  given  a polarization  in  some  plane  instead  of  an 
electric  field.  This  problem  can  be  solved  by  the  following  observa- 
tion. The  relation  between  the  electric  fields  in  two  planes  of  an 
optical  system  is  given  by  a Huygens -Fresne 1 integral.  For  small 
angles,  however,  this  integral  has  exactly  the  same  form  as  the 


1 


39  Uo 

integral  relating  the  vector  potential  to  the  current.  Since, 

spatially,  the  electric  field  is  proportional  to  the  vector  potential 
(the  scalar  potential  is  zero)  and  the  polarization  is  proportional 
to  the  current,  we  conclude  that  the  field  in  the  image  plane  can  be 
found  by  treating  the  polarization  in  the  up-converter  as  an  electric 
field.  Our  method  for  finding  the  electric  field  in  the  image  plane 
will  then  consist  in  first  calculating  the  contributions  to  that  field 
due  to  different  slices  in  the  up-converting  medium,  using  the  results 
of  Fourier  optics,  and  then  adding  these  contributions.  Throughout 
the  calculation,  the  formalism  of  Ref.  39  is  used. 

The  electric  field  in  the  image  plane,  due  to  a slice  of  thickness 
da  at  a distance  Oc  from  the  center  plane  of  the  up-converter,  is 
proportional  to: 


ik  a(x?  + y?) 
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The  quantity  F is  the  transform  of  the  field  in  the  cross  section  a 
of  the  up-converter,  while  A is  the  transform  of  the  Gaussian  aperture 
function.  The  field  in  the  up-converter  can  be  expressed  in  terms  of 
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the  field  behind  the  lens  L, 


i-rrX.  (f  + a)  (xt  + y2) 
1 ' 1 i' 


<x  e 


2 2 

X^F 
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X F, 


'1  \ X F X F / 


(B-3) 


where  F is  the  transform  of  the  field  behind  lens  L . Finally, 
L1  1 
the  field  behind  can  be  expressed  in  terms  of  the  field  in  the 

object  plane: 
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where  M = (X  F/X.f)  . Combination  of  (B.2),  (B.3),  and  (B.M  yields 
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This  result  can  be  greatly  simplified  by  realizing  that  the  phase  fac- 
tor inside  the  brackets  can  only  be  allowed  to  vary  less  than  rr  from 
the  center  of  the  image  to  the  edge,  and  thus  can  be  considered  con- 
stant over  the  aperture  function  A and  taken  out  of  the  brackets. 


I • - 


Hence,  to  an  excellent  approximation: 


(B.6) 


The  total  field  in  the  image  plane  is  found  by  integration  over  Ci 
The  result  is: 
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(B.7) 


We  find  that,  if  there  were  no  change  in  wavelength,  i.e.,  X , 

the  image  would  simply  be  a magnified  replica  of  the  object,  with  a 
loss  in  resolution,  however,  due  to  diffraction  through  the  Gaussian 
aperture.  Because  of  the  change  in  wavelength,  however,  the  image 
gets  modulated  by  a function  that  depends  on  the  length  f the  medium. 
The  effect  of  this  finite  length  will  be  small  if  the  sine  function 
is  close  to  unity  over  the  entire  image.  This  requires: 
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Using 


A 2 

A.  = 7TT. 

1 1 

max 

M = X F/X.F 
s 1 

Aq  - A^M2 

this  can  be  rewritten  as: 


(B.5 


Note  that  the  requirement  BmaxL/2  ;$  n/2  limits  the  change  between 

r =0  and  r.  = r.  of  the  phase  factor  mentioned  above  indeed 

l ii 

max 

to  less  than  n , and  thus  taking  that  phase  factor  out  of  the  con- 
volution was  permitted. 

We  can  now  easily  calculate  the  number  of  resolution  elements  in 
the  image  as  the  ratio  of  the  area  of  the  image  to  the  area  of  a dif- 
fraction spot.  The  area  of  the  image  is  simply  the  area  of  the  object 
times  the  magnification: 


Ai  = Ao^ 

The  area  of  the  diffraction  spot  could  for  example  be  defined  as  the 
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area  between  the  50$  points  of  the  Gaussian,  i.e.. 


x2f2 

a = in  2 t; — 

m,0/2 


where  is  the  area  of  the  pump  beam.  The  resulting  resolution  is: 


A.  1 A A 1 A Ail 

B = - X = -V 

a.  In  2 X7f  In  2 K 

LI  1 


using  f = A . Hence,  approximately: 


Using  the  maximum  value  for  the  solid  angle  of  the  up-converter,  given 
by  (B.8) : 


(B.10) 


The  resolution  found  using  the  above  criterion  is  a factor  of  two 


higher  than  what  is  commonly  used,  i.e.,  R = A > , where  for 

our  case  A = A /2  . This  more  conservative  value  has  been  used  in 

P 

resolution  calculations. 


APPENDIX  C 


CHROMATIC  ABERRATION 


In  Appendix  B we  found  that  the  up-converter  magnifies  the  object 
by  a wavelength  dependent  factor: 

X F 

M = — X - (C.l 

X.  f 

1 

To  find  a criterion  for  chromatic  aberration,  we  will  consider  an  ob- 
ject that  only  has  one  (arbitrary)  spatial  frequency  (in  one  dimension  , 
but  is  illuminated  by  two  different  IR  frequencies.  We  take  the  ampli- 
tudes of  these  two  frequencies  to  be  the  same.  The  field  in  the  object 
plane  can  then  be  written  as 


UQ  = cos  (kxlxQ)  + cos  (kx£?x0) 


with 


xl 


= kx2 
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or 


(C.2) 


in  terms  of  the  angles  0^  and 
k vectors  of  the  two  IR  waves  . 


0,  between  the  z direction  and  the 
The  magnified  field  in  the  image  plane 


is  given  by: 


(c.3 


So  at  the  output  we  have  two  different  spatial  frequencies,  and  there 
will  be  a beating  between  those.  Clearly  what  we  want  is  that  the 
beating  period  is  longer  than  the  size  of  the  image.  Let's  first  re- 
write (C.3)  such  that  the  beating  period  is  shown: 


/ k ,/M,  + k ,/M.  \ / k ,/M,  - k /M  \ 

/ xl  1 x2  2 1 / xl  1 x^  L \ 

I x.  cos  x.  

\ 1 2 / \ 1 2 / 


U.  * cos 
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If  we  limit  the  phase  of  the  beating  factor  to  nf^  , then  the  inten- 
sity in  the  image  plane,  at  the  edge  of  the  image,  will  be  half  the 
value  in  the  center.  So  we  require: 


(kxI/Ml  ' lV/M  'ri  l 


,c  .5 
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If  we  let: 


03.^  = 03.  „ + Zto. 

i2  ll  1 


03  _ = 03  , + &13. 

s 2 si  x 


this  condition  can  be  rewritten  as: 


1 X F 
_ s 

14  fVi 


(C.6; 


for  all  angles  9 < 9^  in  the  object  plane.  Clearly  the  aberration 

max 

is  most  severe  for  large  , i.e.,  large  spatial  frequencies.  Ex- 

pressing the  radius  of  the  image  in  terms  of  the  radius  of  the  object 
and  the  magnification,  the  condition  (C.6)  for  the  maximum  angle  0 
can  be  rewritten  as: 


£03 


03. 


< — 


1 coeo 


(C.7) 


max 


Using: 


= 7TT 


AU0  = 710 

R 


0 

2 

0 

max 
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this  can  be  rewritten  as: 


An  7T  1 

05^  ~ h 


(c.8) 


The  criterion  for  chromatic  aberration  can  then  be  written  as: 


(C.9) 


The  aberration  is  most  visible  near  the  edge  of  the  image,  and  is 
worse  for  higher  spatial  frequencies. 

The  IR  bandwidths,  as  allowed  by  chromatic  aberration,  for  3 P 
and  10  p up-conversion,  and  for  various  resolutions,  are  shown  in 
Table  VII. 
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APPENDIX  D 

PRESSURE  BROADENING  ESTIMATES 


The  formalism  developed  in  Ref.  ^1  was  used  to  estimate  the  pres- 


sure  broadening  of  the  6s-7s  and  7s-7p  transitions  in  cesium.  Matrix 

elements  were  taken  from  Ref.  19. 

A.  TRANSITION  6s" ^ /0  - 7s‘  S1//o 

A B 

Paths  used  for  the  shift  of  6s  + 6s 

Atom  A 

Atom  B 

6s  - 6p 

6s  - 6p 

A 

Paths  used  for  the  shift  of  6s  + 

1 

Atom  A 

Atom  B 

6s  - 6p 

7s  - 6p  1 

6s  - bp 

7s  -7p  I 
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One  finds: 


= 1.1  X 10 


-26  -i 

5 X 10  X N [cm  ] 


where  N is  in  [m  . 


B.  TRANSITION  - 7p2?.^/2 


A B 

Paths  used  for  the  shift  of  6s  + 7s 


Atom  A Atom  B 


6s  - 6p  7s  - 6p 

6s  - 6p  7s  - 7p 


A B 

Paths  used  for  the  shift  of  6s  + 7p 


Atom  A 


Atom  B 


7p  - 6s 


fp  - (s 


6s  - 6p 


7p  - 8s 
7p  - 9s 
7p  - 5d 


7p  - 6d 


7p  - 7d 
7p  - 8d 
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One  finds : 


2 X 10 


-hi 


-2-  = 5-5  x 10"U 

ft 


If  only  the  1 /r'  potentials  are  used,  p,,  = 2b  X . If  the 

(Jm 

and  l/R^  potentials  are  used  with  the  same  sign,  P = 27 

vTTTl 

the  l/Rl  and  l/R^  potentials  are  used  with  opposite  signs, 
21  X . The  resulting  linewidth  varies  between  ^.7  X 10  ^ X N 
and  7.7  X 10  ^ x N [cm  , where  N is  in  [m  ' ]. 


l/R6 
. if 

P0m 

[cm  1 ] 
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APPENDIX  E 


DEPENDENCY  OF  SUSCEPTIBILITIES  ON  POLARIZATION 


A.  INTRODUCTION 

In  this  appendix  we  study  the  effect  of  the  polarization  of  the 
electric  fields  used  in  linear  and  nonlinear  processes  on  their  sus- 
ceptibilities. Although  most  of  the  results  derived  are  intuitively 
quite  obvious,  it  is  interesting  to  see  how  they  are  derived  quantum- 
mechanically . Some  of  the  results  are: 

— The  dipole  moment  set  up  in  a first  order  process  is  parallel 
to  the  exciting  electric  field,  and  has  a magnitude  independent 
of  the  polarization  of  this  field. 

(5) 

— The  magnitudes  of  the  third  order  susceptibilities  ^zzz  anc* 
X(5>  for  a path  s ->  p ~>s  -»  p -»  s are  equal,  i.e.,  it  does 
not  make  a difference  whether  the  electric  field  used  in  the 
third  step  is  parallel  to  or  orthogonal  to  the  field  used  in 
the  first  two  steps.  The  generated  dipole  moment  is  also 
parallel  to  the  field  used  in  the  third  step. 

— If,  in  the  previous  process,  a path  s-»p->d-»p-»s  is  used, 
then  the  susceptibility  for  the  case  of  the  three  fields  parallel 
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is  a factor  of  two  higher  than  for  the  other  case.  This  becomes 
understandable  if  one  looks  at  the  spatial  distributions  of  the 
wave  functions  used  in  both  processes. 

The  expressions  for  the  linear  and  nonlinear  susceptibilities  are 
given  in  Ref.  12.  For  this  discussion  we  are  not  interested  in  consid- 
ering all  the  permutations  of  the  electric  fields,  so  we  drop  the  S 
operator  (=  symmetrization  operator).  Then  we  still  have  to  sum  over 
all  the  atomic  levels.  If  an  atomic  level  is  (2m  +1)  fold  degen- 
erate, it  has  to  be  looked  at  as  consisting  of  (2m  + l)  different 
levels,  all  at  the  same  energy,  and  we  have  to  sum  the  contributions 
of  all  of  them.  Except  for  this  summation,  we  also  have  to  sum  over 
all  the  possible  paths  (e.g..  Is  -» 2p  , but  also  Is  -* 3p  , etc. 
for  a first  order  process  in  H ) . We  will  also  drop  this  last  sum- 
mation because  it  is  irrelevant. 


B.  DEPENDENCE  OF  THE  DIPOLE  MOMENT  GENERATED  IN  A FIRST  ORDER 
PROCESS  ON  THE  POLARIZATION  OF  THE  ELECTRIC  FIELD 


From  Ref.  12: 


X^  (u>) 

pa  ' ' 


(E.l) 
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where 


h 

a 


direction  of  generated  dipole  moment 
polarization  of  electric  field 


(along  x,y,  or  z) 


a,  b 


eR1 


ab 


aa 


atomic  energy  levels 

matrix  element  between  a,b  of  the  p component  of 
the  dipole  moment  operator 

atomic  number  density 

equilibrium  occupancy  of  level  a ; we  set  p°  = 1 

SL  di 

in  what  follows 


For  our  purposes,  as  discussed  in  Section  A: 


(®) 


m-deg.  ab  ab 


(E.2) 


Using  this  expression  for  the  first  order  susceptibility,  we  will: 

(a)  Calculate  the  dipole  moment  set  up  in  the  z-direction  by  an 
electric  field  in  the  z-direction. 

(b)  Calculate  the  dipole  moment  set  up  in  the  x-direction  by  an 
electric  field  in  the  x-direction. 

(c)  Calculate  the  dipole  moment  set  up  in  the  y-direction  by  an 
electric  field  in  the  x-direction. 
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(d)  Prove  that  the  linear  dipole  moment  set  up  by  an  electric 
field  is  parallel  to  that  field,  and  that  its  magnitude  is 
independent  of  the  polarization  of  the  field. 

1 . The  Dipole  Moment  Set  Up  in  the  z-Direction  by  an  Electric  Field 

12 

in  the  z-Direction  is: 


vt; 


/ 


dcoE  (cd)  (cu)  e"1;t 

z zz 


where 


(cd) 

zz  v 7 


■:  ? ( 


, , - 03  03„,  + 03 

m-deg.  01  01 


(E.3) 


from  Section  D: 


(s | z | p,m  = 0)  = i/V? 

(s|z|p,m  = ± 1)  = 0 


so  in  this  case  the  sum  only  contains  one  term;  for  m = 0 : 


m N 1 

XUj(03)  = - - - 

zz  ' 


ft  3 \“01  - “ O>01  + 03 


(E  .ft ) 
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2 . The  Dipole  Moment  Set  Up  in  the  x-Direction  by  an  Electric  Field 
in  the  x-Direction  is: 


CO 

p,(‘)  - / ‘V"’  x»(»>  *'la>t 

-CO 


where 


x(i) 


(a)) 


from  Section  D: 


N 

h 


E 

m-deg . 


|x  |2  |x  |2 
*P . — sp' 


*01  ‘ * *01  + * 


(E.5) 


(s|x|p,m  = 0) 

(s  | x | p,m  = ± l)  = + l/^7P 


hence,  summing  over  m : 


X^(cd)  = 

XX 


- - M(s|x|p,m  = - 1 ) I ‘ + | (s  | x|p,m  = 1>M 


1 1 
+ 


*01  - * *01  + * 


(I+I)(_i L_) 

h \ i 6 / \ui  _ to  , + ai  ! 


01  - D *01 


) (®)  (E.6) 


One  shows  similarly:  X^^  (o>)  = ^ (to)  = x^'(co) 


lOQ  - 


Field 


This  follows  immediately  from  Sections  1,  2,  and  J . A formal 


,(U, 


proof  that  takes  into  account  the  frequency  dependence  of  Xv  ' (o>) 
goes  as  follows: 

Specify  the  electric  field  by  its  direction  cosines,  i.e.. 


— 2 2 2 

i£(E  cos  a , E cos  p , E cos  7)  , with  cos  a + cos  p + cos  7=1; 


then: 


XX 


p (t) 

X 

= cos  a J 

r 

den 

(a>)  E(ffl)  e 

-iu>t 

= COS  P J 

dU5 

X^(cn)  E(cn)  e 

yy 

- i'-Dt 

Pz(t) 

= cos  7 J 

f 

do> 

E (CD)  e 

-iu>t 

^ >> 
X* 
II 

1 (o>)  = X^ 

zz 

(“>) 

, it  follows 

Px(t) 

= 

I COS  CL 

py(t) 

= 

I COS  P 

P2(t) 

= 

I cos  7 
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1 

I 


where 


1(0  = 


f 


dtu  X^1  ^ (oi)  E(cn)  e i nt 


hence 


P(t)//E(t) 


| P ( t ) | = I , independent  of  the  direction 

of  the  electric  field 


(E.9) 


C.  EXAMPLES  OF  THE  DEPENDENCY  OF  THE  THIRD  ORDER  NONLINEAR  SUSCEPTI- 
BILITY ON  THE  POLARIZATION  OF  THE  ELECTRIC  FIELDS 

Since  the  number  of  possible  third  order  processes  is  very  large, 
we  will  only  do  a few  specific  examples.  The  third  order  suscepti- 
bility, as  given  in  Ref.  12,  has  four  terms.  For  convenience  we  will 
only  consider  one  of  them,  more  specifically  the  one  that  lias  three 
minus  signs  in  the  denominator,  which  is  usually  the  dominant  one. 
Again  we  do  not  permute  the  frequencies,  i.e.,  we  apply  the  fields 
in  a well  defined  order,  e.g.:  in  the  first  two  steps  we  apply  a 

field  of  frequency  ni  along  the  z axis,  in  the  third  step  and 
at  frequency  0L  along  the  x axis.  We  also  do  not  sum  over  all 
the  atomic  levels  as  explained  in  Section  A.  So  for  our  purposes 
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(Ref.  12): 


(t) 

M- 


•///*“  l®2'to5X^r(<“l’“2’"3,EC<(<“l>Ef3<"2>S<°>3) 


where 


x(3) 

’ pQp7 


(ai1,a>2,m3 ) 


E 

m-deg . 


Ra  RP  R ' R^ 

ab  dc  cd  da 

(<W  , - 0}  ) (CD  - CD,  - CD  ) (CD  , - CD  - CD  - CD  ) 

v ab  1 ; v ac  1 2 1 x ad  1 2 3 ' 

(E.10) 


Using  this  expression  for  the  third  order  susceptibility,  we  will  do 
the  following  examples: 


For  a path  s -»  ? — > s — > p — > s 

(a)  Dipole  moment  along  z when  all  fields  are  along  z , i.e.. 


(b)  Dipole  moment  along  x when  the  fields  used  in  the  first  two 
steps  are  along  z , and  the  third  field  is  along  x , i.e., 

x<5>  . 

xzzx 

For  a path  s -»  p ->  d -»  p ->  s 

(a)  Dipole  moment  along  z when  all  fields  are  along  z , i.e., 

zzzz 

(b)  Dipole  moment  along  x when  all  fields  are  along  x , i.e., 

x(3) 

xxxx 
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(c)  Dipole  moment  along  x when  the  fields  used  in  the  first  two 


steps  are  along  z , and  the  third  field  is  along  x , i.e., 

x(5> 

xzzx 

To  shorten  the  algebra,  we  will  represent  the  resonance  denominator 

(“oi  ~ V (m02  • mi  • “2}  (6D03  - ®1  - ®2 " “3 } by  one  say  D • 

1 . Process  Using  a Path  s->p-»s~>p->s 
(a)  Calculation  of  X^5  ' 


- y 


z z z z 


m-deg . 


From  Section  D: 


(s|z|p,m  = 0)  = 

(s|z|p,m  = ± l)  = 0 


hence,  summing  over  m : 


N / i \ " 1 N / 1 W 1 \ 

6ft5  \ \fT  I D 6ft5  \ 9 / \ D / 


(b)  Calculation  of  X' 


- -X  y 

xzzx.  t—j 


Z Z X X 


m-deg . 


(E.13) 
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r 


From  Section  D: 


(s|z|p,m  = 0)  = l/^T 

(s|z|p,m  = ± l)  = 0 


(s | x ] p,m  = 0)  = 0 

(s|x|p,m  = ± 1)  = + l/^o"1 


hence,  summing  over  m : 


X 


(3) 

xzzx 


N 1 

6ft5  D 


(s|z|p,m  = 0) {p,m  = o|z|s) 

(s|x|p,m  = - l)(p,m  = - l|x|s) 


X 


+ (s|x|p,m  = l)(p,m  = l|x|s) 


N 1 

6 tP  D 


6 6 


N / 1 \ / 1 \ 

\ q / V d / 


,(3) 

zzzz 


( E . 1 V 


So  it  does  not  make  a difference  in  this  case  whether  the  third  field 
is  parallel  to  or  orthogonal  to  the  first  two.  This  could  be  expected 
since  the  vibration  set  up  in  the  second  step  (s  level!  is  spatially 

/ 2 \ / ^ \ 

symmetric.  It  can  also  easily  be  shown  that  X ' and  X1,  are 

zzzx  yzzx 

zero,  i.e.,  the  generated  dipole  moment  is  parallel  to  the  third 
field . 
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k 


2. 


Process  Using  a Path  s p ->d  -»  p «->  s 


(a)  Calculation  of 


x(5) 

zzzz 


I 

l 


x(3) 

zzzz 


£ 

m-deg. 


z z ,z,  z 
sp  pd  dp  ps 

D 


(E.15) 


N / 1 


xxxx  6ft5  \d 


(s|x|p,m=-l) (p,m=-l|x|d,m=0) (d,m=0 | x| p,ra^-l) (p,m=-l |x|s) 

+ (s | x|p,m=-l) (p,m=-l |x|d,m=0) (d,m=0 | x| p,m=l> (p,m=l |x|s) 

+ (s | x| p,m=-l) (p,m=-l|x| d,m— 2) <d,m=-2|x|p,m=-l) (p,m=-l |x|s) 
+ (s |x|p,m=l) (p,m=l|x|d,m=0) (d,m=o|x| p,m=-l> (p,m=-l |x| s) 

+ (s | x | p,m-l) (p,m=l |x|d,m=0) (d,m=o| x| p,m=l) (p,m=l |x| s) 


(s|x|p,m=l) (p,m=l|x|d,m=2) (d,m=2 | x| p,m=l) (p,m=l | x| s) 


Using  the  values  for  the  matrix  elements  given  in  Section  D: 


(E .18) 


N / 1 


GfP  \ D / 

( i ) ( - 1 ) ( - 1 ) ( i /^p ) 

+ d/vr)(-i/V5o")(1/V5°")(-i/>F) 

+ (i/^P)(i/^T’)(i/^51)(i/^61) 

+ (-i/VT  ) ( ) ( 1/^jo1 ) ( - 1 /^r ) 


(Cont'd) 
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The  susceptibility  for  the  case  where  the  three  fields  are  parallel  is 
a factor  of  two  higher  than  in  the  case  where  the  third  field  is  or- 
thogonal to  the  first  two.  This  is  not  too  surprising  if  one  looks  at 
the  spatial  form  of  the  wave  functions  involved. 

D.  MATRIX  ELEMENTS  OF  THE  DIPOLE -MOMENT  OPERATOR 

Since  only  first  and  third  order  processes  are  considered,  only 
s , p , and  d levels  are  involved,  i.e.,  i = 0,1,2  . For  an 

electron  moving  in  a spherically  symmetrical  potential  the  wave  func- 
tions can  be  written  as  (Ref.  hj>) : 

( l = 0,1,...,  n-1 

*(r,0,cp)  = R ,(r)Y">,cp)  J (E  .23) 

( m = -l  , -f+1  , . . . , f-1  , l 

The  Y^(0,q>)  for  1 = 0,1,2  are  (Ref.  b2) : 


(3  cos'  0-1)  (Cont'd) 
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sin  9 cos  9 


1 e ) 


The  matrix  elements  of  the  position  operator  are  defined  by: 


(nim|~r|n'i'm') 


where  |nfm)  = ,(r)Y™(0,cp)  . Since  for  this  discussion  the  value 
of  the  radial  integral  is  irrelevant,  we  will  set  its  value  equal  to 
one.  So  we  have  to  evaluate  the  following  integrals: 


( £m| cos  0 | I'm') 

(£m|sin  0 cos  cp|i'm') 

(£m|sin  0 sin  Cp|£'m') 

The  Ym(0,cp)  are  given  by: 

Y^(0,cp)  = P™(cos  0)  e1”1^ 

From  this  follows  the  selection  rule  for  £ : Af  = ± 1 . Since  for 

this  discussion  the  selection  rule  for  m is  important,  we  will  derive 


for: 
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it  here.  The  cp  part  of  the  matrix  elements  is  of  the  form: 


( e 


imC,’ 


27T 

eim>  ) = f 

0 


-imcp  / _ I im'cp 

e I cos  cp  I e dcp 

, sin  cp  ; 


hence : 


2n 


f < 


icp  -icp. 
e + e ) 


/2 


im'cp 


- e-lCP)/2 


c(e 


for: 


;.^n  = 0 


■ \n  = ± 1 


y : :^m  = ± 1 


List  of  relevant  matrix  elements: 


(l)  s p 

(s  I z|  p,m  = o)  = 

(s | z | p,m  = ± l)  = 0 

(s | x | p,m  = 0)  =0 

(s  | x | p,m  = ± 1)  = + l/^i”1 


dcp 
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r 


/ 

vs 

|ylp,m  = 0) 

= 0 

/ 

\S 

|y|p,m  = ± 1)  - - 

■ i/-^* 

2_->_d 

(only  the 

non-zero  ones  are  given) 

(p,m 

= 0 1 z ] d,m  = 

0)  = 

(p,m 

= ± 1 | 2 |d,m 

= ± 1) 

- i/vr 

(p,m 

= o|x|d,m  = 

± 1)  = 

= + i/^To' 

xp,m 

= - 1 |x| d,m 

= 0)  = 

= - 

(p,m 

= - 1 1 x| d,m 

= - 2) 

= i/^T 

(p,m 

= 1 j x ( d,m  = 

0>  = 

l/yJW 

(p,m 

= 1 | x| d,m  = 

2)  = 

- l/yfT 

\P,m 

= 0|y|d,m  = 

± 1) 

= ± i/^To1 

\P>m 

= - l|y|d,m 

= o> 

- i/yprt' 

(p,m 

= - 1 | y | d,m 

= - 2) 

= i/^T 

\ p,m 

= l|y|d,in  = 

0)  = 

- i/yfW 

\p,m 

= l|y|d,m  = 

2>  = 

i/^T 
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APPENDIX  F 


CONVERSION  OF  LINE  STRENGTHS,  OSCILLATOR  STRENGTHS 
AND  EINSTEIN  COEFFICIENTS  TO  MATRIX  ELEMENTS 


USED  IN  PERTURBATION  CALCULATIONS 


A . INTRODUCTION 


In  this  appendix  we  address  the  problem  of  calculating  the  matrix 
elements  to  be  used  in  each  step  of  a perturbation  sequence,  given  the 
relevant  line  strengths,  oscillator  strengths  or  Einstein  coefficients. 
Since  in  practice  we  are  almost  always  concerned  with  transitions  s-p 
or  p-d  , only  these  are  discussed.  We  consider  two  cases: 

Case  1 : The  frequency  detuning  is  much  larger  than  the  term  split- 

tings of  the  multiplet.  In  this  case  one  can  use  the  to- 
tal matrix  element  for  the  transition  (n l)  - (n'f' 

Case  2:  The  frequency  detuning  is  comparable  to  or  smaller  than 

the  term  splittings  of  the  multiplet.  In  this  case  one 
has  to  take  into  account  the  fine  structure  of  the  multi- 
plet, i.e.,  one  has  to  calculate  the  matrix  elements  for 
each  of  the  transitions  (n£j)  - (n'f'j')  , and  after- 

wards sum  over  j and  j ' . The  results  in  this  appendix 
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also  show  how  to  distribute  the  total  line  strength  of  a 


multiplet  among  its  various  lines. 

Clearly,  if  one  neglects  the  energy  difference  between  the  levels  of 
a term,  after  summing  over  j and  j'  in  the  second  case  one  should 
find  the  same  final  result  as  in  the  first  case.  We  demonstrate  this 
by  calculating  a third  order  nonlinearity  in  both  schemes,  after  first 
considering  the  effect  of  spin. 

All  the  results  are  valid  for  atoms  with  ne  electron  in  the  outer 
shell,  for  which  the  unperturbed  wave  functions  are  products  of  a ra- 
dial part  and  a spherical  harmonic.  Also,  the  electric  fields  used  in 
the  perturbation  sequence  must  all  have  the  same  polarization;  if  this 
is  not  the  case,  the  combination  of  the  results  of  Appendix  E and  the 
techniques  illustrated  here  will  still  lead  to  the  answer,  but  each 
case  has  to  be  solved  separately. 

B.  CALCULATION  OF  MATRIX  ELEMENTS  FROM  TOTAL  LINE  STRENGTHS, 
DISREGARDING  TERM  SPLITTINGS 

Remarks : 

- The  line  strength  S is  defined  between  two  levels  (Sections 
G and  H),  e.g.,  (s,j  = l/2)  - (p,j  = j/2)  . However,  if  we  sum  the 

line  strengths  of  the  transitions  (nij)  - (n'l' j')  over  all  possible 
pairs  (j,j')  , then  this  total  S is  by  definition  z|(aiF|b)|^ 

between  the  configurations  (nf)  and  (n'f')  (usually  it  is  this 
total  S which  is  given  in  tables).  From  this  total  line  strength 


- 211*  - 


we  will  find  the  matrix  element  to  be  used  in  the  transition  (nf)  - 
(n '/')  . 


— The  eigenstates  used  in  this  calculation  are  n£sm  m ) states, 

i s 

Since  s - l/2  and  n is  irrelevant  for  this  discussion,  we  will 
usually  omit  these  two  quantum  numbers.  The  numerical  values  used  in 
the  following  calculations  are  taken  from  Ref.  U5 . 

— Selection  rules: 


± l 


1 . s 


Transitions 


1 -1 


1 -1 


m =l/2 

s 


m = - 1 
s 


S = e^ 


E (i- 


ly.J  + i2,' 


using  2sn^  = 0 : 


I 


2 

— Calculation  of  L|zab|^ 

= 2|(s  0 l/2 | z I p 0 1/2)  1 2 = 2 (1/3  : = 2 h 

2 

— Calculation  of  Zlx^T  using  ~*n  . = ± 1 : 

2=  2 [ I ; s 0 l/2 1 x[  p 1 l/2)  | + I (s  0 l/ 2 1 x | p - 1 l/2>!2] 

= 2 (l/6  +l/6)  = 2/3 

hence 

Z)ixabr  = £!yaJ2  = S^Zab!2 

as  it  should  be.  As  a consequence  S must  be  divided  by  3 if  lin- 
early polarized  light  is  used.  We  will  see  in  Section  D that  another 

factor  of  2 appears  because  of  the  equal  distribution  of  the  atoms  in 
the  first  level  used  in  the  process  between  the  two  spin  states. 
Example:  Lithium 


2. 


E. 


d Transitions 


1-12-2 


1-12-2 


1 -1 


0 1 -1 


m = 1/2 

s 


m = - l/2 
s 


Calculation  of  ^lza^  using  = 0 : 


7>  2 | I'p  0 l/2 1 2 1 d 0 l/l  I2  + I pi  1/2  1 z I d 1 1/2)1“  J 
| + |(p  - 1 l/2 | z | d - 1 1/2) | - j 

= 2 (U/15  + 1/5  + 1/5)  = 2 (U/15  + 6/15)  = U/j 


Calculation  of  Z|xa^|  using  = ± 1 : 


|<P  0 l/ 2 | x | d 1 1/2) I2  + I(p  0 1 /2 1 x I d - 1 l/2>|2 

^ = 2 + I (P  1 l/2|  x|  d 0 l/2)  |2  + I (p  1 l/2 1 x|  d 2 l/2)  |2 

|{p  - 1 l/2 | x| d 0 l/2) | 2 + |<p  - 1 l/2|x|d  - 2 l/2) | 2 

= 2 {(1/10  + 1/10)  + (1/30  + l/5^  + (1/30  + 1/5)} 

= V3 
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hence 

2>j2  - 5>j*  £'*J 

and  S must  be  divided  by  3 f°r  linearly  polarized  light. 

Now  consider  a nonlinear  process  starting  in  a s-level,  and  with 
a transition  p -> d in  some  intermediate  seep,  e.g.,  in  the  second 
step;  assume  also  that  the  electric  field  is  polarized  along  the  z 
axis  (for  convenience): 


0< 

> ll 

1 

1 

1 

-li 

1 

2 

oJ 

1 1 

-1 

I 

0 

m =l/2 


Identical 

Diagram 


m - 1.  . 


Because  of  the  selection  rule  , not  all  the  's  can 

be  used  in  going  fnm  p to  d The  fraction  used  is: 

(°->°)  _ Vl3 

(0  ->  0)  + (1  -»  1)  + -1  "I  U/l5  + t /l‘ 
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As  a consequence,  in  a perturbation  sequence  step  as  this  one,  the 
line  strength  S has  to  be  divided  by  .5,  besides  the  factor  of  / 
because  of  polarization.  Also  again  an  extra  factor  of  appears 
because  of  spin  as  will  be  explained  later. 

Example : Lithium 


Transition  1 

s 

Vs/15 

z Ref.  19 

Difference  in  j 

2p  - 3d 
1 2p  - Ud 
5p  - 5d 

. • 

11.1 

I 

. ■ 

. 

• 

1.U  j 

, 

1.-5 

• 5 
1. 

C.  CALCULATION  OF  MATRIX  ELEMENTS  BETWEEN  PAIRS  OF  LEVELS  OF  TWO 
TERMS,  GIVEN  THE  TOTAL  LINE  STRENGTH  OF  THE  MULTIPLET 

Remarks : 

— From  a mathematical  point  of  view,  one  can  always  transform 
from  one  complete  basis  set  to  another  one,  in  our  case  from  n . sm  .m 
to  nisjm.  . Clearly  the  magnitude  of  physical  quantities  calcu- 
lated in  both  schemes  should  be  the  same.  The  reason  why  the  n.sim 

J 

scheme  is  a good  scheme  to  use  in  the  presence  of  spin-orbit  interac- 
tion is  that  the  total  Hamiltonian  H + H.  _ is  to  an  excellent  ap- 
proximation  diagonal  in  this  scheme  Ref.  •••'•  . As  a consequence  the 
L-S  interaction  does  not  mix  states  with  different  j values,  but 
only  gives  them  different  energies.  If  this  had  not  been  the  case. 


lq 


the  matrix  elements  calculated  between  two  pure  j -levels  would  have 
been  incorrect  for  real  atoms  where  the  L-S  interaction  is  impor- 


— To  calculate  matrix  elements  between  , n<sjnK  eigenstates, 
one  has  to  express  them  in  terms  of  (nfsm.m^  eigenstates.  The 
coefficients  of  this  linear  transformation  are  the  Clebsch-Gordan 
coefficients.  They  are  tabulated  in  Section  1 for  the  cases  of  in- 
terest to  us,  i.e.: 


= l/2 


0,1, 


— A useful  property  of  matrix  elements  in  the  jnu  scheme, 
that  follows  immediately  from  a symmetry  argument,  is: 


Ij111  j M f'j  "nij  ' ! 


| (l j - m. ] r | i ' j ' - m' ' 
J J 


We  will  frequently  make  use  of  it  to  simplify  the  calculations. 

— We  will  use  spectroscopic  notation,  but  omit  the  superscript 
s -t-  1 2 for  convenience,  i.e.,  we  will  write,  e.g.,  P..  , in- 


stead of  P 


3/2  • 
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— The  selection  rules  in  this  scheme  are: 


Since 


1 . s 


!Aj  = 0 , ± 1 (0  — 0 forbidden) 

An.  = 0 , ± 1 

m.  = m,  + m , An.  = An.  + An  = An 
j l s j l s £ 


Numerical  values  for  matrix  elements  are  taken  from  Ref 


- p Transitions 


Levels 

Spectroscopic  Notation 

s j = l/2 

2c 

Sl/2 

j - 1/2 

Pl/2 

P 

. • /o 

2„ 

J - 3/2 

* 1 

The  transitions  allowed  bv  the  selection  rule  on  i are: 


Sl/2  ' Pl/2 


Transitions 


a . 


Pl/2  l/2  -1/2 

Sl/2  l/2  ^172" 

. {2 

— Calculation  of  L z , using  An.  = 0 , An  =0 

ab  0 j 7 s 

y*i  - 2 i(s  1/2 1/2 1 z 1 P 1/2 1/2)  p 

= 2 I - ^1/3'  (s  0 l/2|z|p  0 1/2) 

+ -^273"  (s  0 l/2 1 z | p 1 - 1/2)  |; 

= 2/9 

The  second  term  in  the  sum  vanishes  because  of  An  - 

s 

— Calculation  of  ^‘lxaiJ  using  - ± 1 , \n 

T.  = 2 | (s  l/2  l/2 1 x|  p 1/2  - 1/2)  |2 

2 I - ^2/3  ‘ (s  0 l/2 1 x | p - 1 1/.  > |C 
= 2/9 
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hence 

2>J2  *2>J*  -E'^l2 

as  in  the  other  scheme,  and  of  course  as  it  should  be.  As  a conse- 
quence S must  be  divided  by  3 for  linearly  polarized  light. 


Sl/2  1/2  -1/2" 

p 

— Calculate  “lzaj,|i'  using  2!m.  = ) , 2tng  - 0 : 

Y-  2 |(s  1/2  1/2 1 z | P 3/2  1/2) |2 

2 I ^2/3'  (s  0 l/2 1 z | p 0 1/2)  I2 

One  shows  again 

£l«J2  E’ZaJ 

and  hence  S must  be  divided  by  3 for  linearly  polarized  light. 
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Conclusion 


Comparing  the  results  of  (l)and(2),  we  see  that  the  line  strength 
for  the  • transition  is  twice  as  large  as  the  one  for  the 

• - P.  / transition.  Hence  the  total  line  strength  is  divided 

among  the  two  lines  of  the  doublet  as  follows: 


2. 


d Transitions 


2 

— Calculation  of  Llz^l  using  /!^n.  = 0 , •/^ns  = 0 

y 2 !(p  1/2  l/2|z|d  5/2  1/2) I2 

= 2 | ^2/5'  ^l/5'  (p  0 l/2 1 z | d 0 1/2) 

+ ^3/V  (P  1 - l/2|z|d  1 - l/2>  i ‘ 

= 2 (26/225  + 2V225)  = V9 

One  shows  again  that  S must  be  divided  by  5 for  linearly  polarized 
light . 

Consider  a nonlinear  process  starting  in  a S^/„  level,  and  with 
a transition  - - D /0  as  some  intermediate  step.  Assume  that  the 

electric  field  is  along  the  z axis: 


°5/2 


Pl/2 


Sl/2 


. 1/2 

1 l/'c 

372 


-1/2 


-1 


-1  (2 


-3 


o 

cl 


In  this  case  of  a transition  - D.  / , the  line  strength  must  not 

be  divided  by  an  additional  factor  because  of  the  fact  that  to  = 

j 

for  all  the  steps  of  the  sequence.  This  will  not  be  so  for  the  transi- 
tions P3/2  - D5/2  and  P?/2  - D5/2  . 


22c 


b. 


P.  /„  - D-,  / Transitions 
i/£ ILL — 


Dj/2  1/2  3/2  -l/2  -3  A 


P3/2  1/2  3/2  -1/2  -3/2 


— Calculation  of  . z | using  2sn.  = 0 , 2*n  = . 

1 ab1  j s 


2 


(|(p  3/2  1 /2 | z 1 d 3/2  1/2) |2  + | p 

( I * ^2/5  ^2/3  p 1 1 z J d 

< + ^3/T  JUT  P 1 - l/2 1 z | d 

+ I - (p  1 1 /2 1 z ! d 1 1 ' 

(1/225  + 9/225)  = ^A5 


3/2  3/2 1 z | d 3/2  3/2)  | 


Again  a factor  of  3 comes  in  for  linearly  polarized  light. 

Consider  a nonlinear  process  starting  in  a S ^ j level,  and  with 

a transition  P / - D . as  some  intermediate  step.  Assume  that  the 
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Calculation  of  fz  , !'  using  vi . , 'ri 

1 ab 1 j s 


E“  2 ' 


p • l/2|z|d  ‘ 1 


V V 


1 , izld  1 


|z|d  5/2  3, 


^ ^ L p 1 - 1 >,11-1 

1 J P 1 1 ' j z ' d 1 1 / 


;51*/225  + 30/225  -./5 


Again  a factor  of  comes  in  for  linearly  polarized  light. 

Consider  a nonlinear  process  starting  in  a level,  and  with 
a transition  P - D as  some  intermediate  step.  Assume  that  tin 
electric  field  is  along  the  z axis: 


1 1 


The  fraction  used  in  this  case  is: 


i-1  *- 


1/::  * 1/ 


hence  the  total  line  strength  of  the  multiplet  is  divided  among  the 
three  lines  as  follows: 


Line 

1 - 

Fraction  j 

Pl/:  -D;/: 

„ i 

P3/:  ' V: 

1/15  1 

p /r  - V; 

9/15 

This  result  agrees  with  numerical  examples  taken  from  Ref. 
Examp le : Na  p-  d 


Trans i t ion 

Line  Strength 

Fraction 

P‘d 

i 

V'  - D 

5/15 

P /.  - D 

5/c  5/' 

5/lr’ 

P - D 

l . 

0/15 

D.  EFFECT  OF  SPIN  ON  THE  CALCULATION  OF  MATRIX  ELEMENTS 
FROM  LINE  STRENGTHS 

As  was  mentioned  in  the  Introduction,  we  require  that  all  the 
electric  fields  used  in  the  process  have  the  same  polarization;  for 
mathematical  convenience  we  take  them  all  along  the  z axis  of 
course  this  does  not  have  an  effect  on  the  final  result  . The  dia- 
grams in  Sections  B and  C may  help  in  clarifying  the  discussion  in 
this  section. 


1.  (ta.m  ) Scheme 
( s 


In  the  calculation  of  a linear  or  nonlinear  process,  one  alwavs 

has  to  perform  a sum  over  all  possible  paths.  In  the  im  ,m  scheme, 

this  sum  can  be  subdivided  into  two  parts:  one  part  with  spin  up 

m = 1/2  and  one  with  spin  down  m - 1 . Because  of  the 

s r s 


selection  rule  'm 


these  two  parts  are  uncoupled.  The  expres- 


sions for  susceptibilities  and  related  quantities  can  then  be  written 


as  : 


( 


fraction  of  the  population  with  ti 
in  the  first  level  of  the  process 


l/!) 


El produc t ' f mat  r ix  \ 
lelemcnts  I 

paths 


Cent  1 d 


(fraction  of  the  population  with  id  - 1 j 
in  the  first  level  of  the  process  ' / 

Z /product  of  matrix! 

lelements  ) 


paths 


Since  the  electronic  matrix  elements  in  both  terms  of  this  sum  are 


equal,  we  can  write  equivalently: 


(total  population  in  the  first  level  of  1 
the  process  / 

E /product  of  matrix! 

lelements  / 


paths 

This  means  that  equivalently  we  can  assign  all  the  population  of  the 
first  level  to  either  one  of  the  spin  states,  and  consider  the  other 
one  empty,  i.e.,  we  only  use  one-half  of  the  diagrams  in  Section  1. 
Because  the  line  strength  includes  a summation  over  spin,  since  it  is 
defined  in  the  (fjm  ) scheme,  this  implies  that  for  each  step  of  a 
perturbation  sequence  the  line  strength  must  be  divided  by  two  to  find 
the  effective  matrix  element  for  that  step. 


( i jm  Scheme 


In  this  scheme  we  cannot  in  general  make  a subdivisi  i as  in  thi 

other  sch< me,  i.e.,  stat  1 • tl  tlues  a 


st  uncoupled,  e.g.,  s 1 1 |x|p  l 


. However,  if  w< 


assume  all  the  electric  fields  polarize  1 al  n>  the  • i%i-  . then  stat< 


with  m.  values  of  opposite  s i -m  art  uncoupled  tn . 


and  wc 


'I 


can  make  an  argument  analogous  to  the  one  used  in  the  other  scheme,  i.e.: 

^ / fraction  of  the  population  with  total  angular 

— l momentum  up  in  the  first  step  of  the  process 

(product  of  matrix 
elements 

paths 

(fraction  of  the  population  with  total  angular 
momentum  down  in  the  first  step  of  the  process 

(product  of  matrix 
elements 

paths 

The  rest  of  the  argument  is  analogous  to  the  previous  case,  and  the  con- 
clusion is  the  same. 

E.  SUMMARY  OF  RULES  FOR  CONVERSION  FROM  LINE  STRENGTHS  TO 
MATRIX  ELEMENTS 

Remarks : 

— Only  the  absolute  value  of  the  matrix  elements  can  be  found  from 
the  line  strength;  in  the  tables  below  the  conversion  is  given  from  line 
strength  to  the  square  of  the  matrix  element. 

— For  p-d  transitions,  we  have  to  consider  ,two  cases  see  Sec- 
tions B and  C for  an  explanation  : 

Case  1 : None  of  the  preceding  steps  in  the  perturbation  sequence 

uses  a s-p  transition. 

Case  . : At  least  one  of  the  preceding  steps  uses  a s-p  transi- 


- 


■ if 


t ion . 


CALCULATION  OF  A THIRD  ORDER  SUSCEPTIBILITY  IN  THE 


F . 

■ im  ,m  AND  ijm.  SCHEMES 
is  j 

The  purpose  of  this  section  is  to  demonstrate  with  an  example  that, 
if  j-degeneracy  is  assumed,  susceptibilities  calculated  in  both  schemes 
have  the  same  value,  as  they  should.  One  must  assume  j-degeneracy  to 
make  the  detunings  for  the  various  transitions  between  levels  of  two 
terms  identical.  From  a practical  point  of  view  it  follows  that,  in 
the  case  of  detunings  much  larger  than  the  doublet  splittings,  suscep- 
tibilities can  be  calculated  in  either  scheme,  the  calculations  in  the 

in  m scheme  being  much  simpler  however.  This  calculation  is  also 
i s 

a numerical  check  on  the  results  of  the  previous  sections. 

-Path:  s-p-d-p-s 

— We  assume  an  electric  field  parallel  with  the  z axis,  hence 

' - y!  s | z | p ^ < p | z | d ' ’ d | z | p ) p j z | s ' 

paths 

— Symbols : 

total  line  strength  for  the  s-p  multiplet 


S 


total  line  strength  for  the  p-d  multiplet 


\ 

Is 

VV 

S^.,/10 

Vs! 

> 

Is  1 \ 

FJi" 

s1s2/i6eo 

Vsi  Vs 

> 

lS  - ' 

SjS  ./l 50 

x a SjSg  1/3-  **  + 

1 1- 

• 1 1 

+ l/lt ; 

+ 1/1501 

= s s ),!  i 

sisr/> 

uu 

DEFINITIONS 

— Level : Set  of  .'j 

+ 1) 

states , 

each  characterized  by 

ferent  m 

value 

. If  the 

atom  is 

unperturbed 

ternal  fields,  the  (;'j  + 1 states  are  degenerate. 

When  the  atom  is  perturbed  they  can  have  different 
energies . 

— Line:  Radiation  associated  with  all  possible  transitions  be- 

tween the  states  belonging  to  two  levels.  The  radiation 
resulting  from  a transition  between  a particular  pair  ol 
states  is  called  a compone nt  of  the  line. 

—Configuration:  The  quantum  numbers  n,  t are  said  to  specify 

the  configuration  of  the  electron. 
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— Term:  The  (one  or)  two  levels  into  which  each  configuration  is 

split  because  of  spin-orbit  interaction  are  said  to  con- 
stitute a doublet  term. 

— Multiplet:  Ensemble  of  all  lines  connected  with  transitions  from 

one  term  characterized  by  L,S  to  another  characterized 
by  L',S'  . 

H.  LINE  STRENGTH,  OSCILLATOR  STRENGTH  AND  EINSTEIN 

Ik 

COEFFICIENT 

I . Line  Strength 

Given  two  levels  A and  B , each  consisting  of  several  states 

( respectively  i.  + 1 and  2j_  +1  ),  the  line  strength  is  defined 
r A D 

as : 

S (A,  B = S B,  A S ! a [ P | b |' 

a y b 

where  P = - e £ r.  is  the  dipole  moment  operator.  If  only  one  elec- 
i 1 

tron  is  involved,  P = - er  and 

S = e2  £ CUJ;?  * IvJ'  1 hj2> 

a,  b 

The  summation  is  over  all  states  a,b  of  the  levels  A.B 


Example : 


P>j  = 3/2 


372"  1/2"  -“1/2  -3/2 


<-  ■ states 


s,j  = l/2 


TTi? 


<—  states 


Note  that  S is  symmetrical  to  the  initial  and  final  levels.  The 
physical  origin  of  this  definition  of  S is  that  the  total  intensity 
of  a line  in  natural  excitation  equal  numbers  of  atoms  in  all  the 
states  of  the  same  level  is  given  by: 


I (A,  B 


t !-7|  v" 

N a hv  r 

j?hc 


E 


I a!p!b)|2 


~ N (a  S( A,B 


Oscillator  Strength 

Given  two  levels  ^aj  and  (a'j'  , the  oscillator  strength 

is  defined  as: 

m v(a'j',aj  S,uj,u'j' 
e h 


f t'  j j 


2j  + 1 


where 


(Cij  ' = lower  level 

a'  j ' ) = higher  level 

v = frequency  of  the  transition 

Note  that  the  oscillator  strength  is  defined  between  two  levels,  each 
having  a particular  j value.  The  physical  origin  of  this  definition 
of  oscillator  strength  is:  it  is  the  ratio  of  the  contribution  of  the 

spectral  line  to  the  atomic  polarizability,  and  the  polar- 

izability of  a classical  harmonic  oscillator  with  a resonant  frequency 
equal  to  the  frequency  v(a'j',aj'  of  the  transition.  Using  more  com- 
mon notation: 

(aj  = i 
(a'j')  = k 

2J  + 1 = 8i 

v(a'j',aj ) c/x 
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and  substituting  numerical  values  for  the  constants: 

(F.l) 

(F.2) 

where 

- X in  X 

S in  a.u.  (i.e.,  in  units  of  a e 
conversion  to  Int . System: 

M = 8.b6  x 1 o'5°  ^1T 

Expressions  (F.l)  and  (F.2)  agree  with  those  given  in  the  conversion 
table  in  Ref.  X5 . Note  that  the  oscillator  strength  is  not  symmetri- 
cal in  the  initial  and  final  levels,  since  S is  symmetrical  and  g^ 
refers  to  the  lower  level.  However,  if  f^  is  known,  one  can  always 
find  S using  (F.2).  As  a consequence  it  was  sufficient  to  give  all 
the  results  in  Section  E in  terms  of  S 

3 . Einstein  Coefficient 

Given  two  levels  A and  B , the  Einstein  coefficient  is  de- 
fined as: 


- 2ltl 


G(A,B) 


uUtt* \r>  S (A,  1 


3hc^  2j  A + 1 


A = initial  level  = higher  level 


B = final  level  - lower  level 


Just  as  the  oscillator  strength,  an  Einstein  coefficient  is  defined 


between  two  levels,  each  having  a particular  j value.  The  physi- 


cal origin  of  this  definition  of  Einstein  coefficient  is:  the  total 


intensity  of  the  line  A — > B is  given  by: 


I A, B = N A hv  ( A, B 


where  N (A  is  the  total  population  of  level  A the  reader  is  urged 


to  read  pp.  and  in  Ref.  at  this  point  to  get  the  complete  pic- 


ture . Using  more  common  notation: 


A - k 


B = i 


jA  ' 1 8k 


Ik 


r 


and  substituting  numerical  values  for  the  constants: 


where 

- i n X 

S in  a .u . 


■ ss Lons  F.;  and  F J agree  with  those  given  in  the  conversion 
table  in  Ref.  5.  Just  as  in  the  case  of  the  oscillator  strength,  the 
Einstein  coefficient  is  not  symmetric  in  the  initial  and  final  levels. 
However,  given  G , one  finds  S using  F.  , and  again  there  is  no 
loss  in  generality  by  expressing  everything  in  terms  of  S 

I.  CLEBSCH-GORDAN  COEFFICIENTS 

One  goes  from  the  basis  set  In  Asm  m to  the  basis  set  jnfsjmj 
by  a linear  transformation.  The  coefficients  used  in  this  transforma- 
tion are  th<  ch-Gordai  ef  1 ts . Cleai  th<  bei  E i en- 

states  in  the  second  scheme  has  to  be  equal  to  the  number  of  eigenstates 
in  the  first  scheme,  since  all  the  eigenstates  are  orthogonal  in  the 
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first  scheme;  this  can  be  shown  very  easily; 


V 


!m,m  basis:  degeneracy 

' i s 


| jrrr  basis:  degeneracy 


2 X (21  + 1 

spin  m degeneracy 

E < i 

m.  degeneracy 
J J 

2.(1  - l/  ■ 1 • 1 ■ l 


2(2 1 + 1 


The  coefficients  given  in  the  following  table  were  calculated  using 
Table  1 on  p,  i'o  in  Ref.  U* . / 

Transformation  Matrices 
1.  s-state 


l = 0 


nl  = 


i 1 ' 


/ ' 1/  1/2  \ /l 

V 1 1/2  -!/:■)/  \ 


m 

s 


1 


-1 


■ 1 


\ 
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